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1.  INTRODUCTION: 


Prostate  cancer  is  heterogeneous,  with  diverse  outcomes.  It  is  unclear  to  what  extent  low  grade, 
indolenct  cancer  can  evolve  into  aggressive  disease.  And  in  aggressive  disease,  it  is  unclear  if  the 
genetic  alterations  more  common  in  late  disease  are  present  early  on,  but  at  low  frequency,  or  if  they 
only  develop  late  in  disease  progression.  This  project  aims  to  characterize  the  evolution  of  both 
localized  and  metastatic  prostate  cancer.  One  aim  seeks  to  evaluate  the  evolution  of  localized  prostate 
cancer  to  high  grade  or  aggressive  variants.  The  second  aim  seeks  to  characterize  the  evolution  from 
localized  to  metastatic  prostate  cancer. 

2  KEYWORDS: 

Cancer  genetics,  tumor  evolution,  tumor  heterogeneity,  prostate  cancer,  exome  sequencing 


3.  OVERALL  PROJECT  SUMMARY  BY  SPECIFIC  AIM: 

The  research  project  proposed  in  this  award  was  comprised  of  two  specific  aims,  with  further  division 
into  sub  aims,  and  tasks  outlined  in  the  statement  of  work.  This  portion  of  the  report  will  give  a  brief 
overall  summary,  followed  by  specific  reference  to  the  tasks  outlined  in  the  SOW. 

Aim  1:  Investigate  the  molecular  relationship  between  coincident  prostate  cancer  foci. 

We  set  out  to  characterize  the  molecular  relationship  between  coincident  foci  in  two  scenarios: 
coincident  low  and  high  grade  disease;  and  coincident  acinar  and 
ductal  histologies.  In  the  first  year  we  published  that  coincident 
high  and  low  grade  prostate  cancer  diverge  early  in  their  evolution, 
with  only  a  minority  of  mutations  being  clonal  (VanderWeele  et  al. 

Cancer  Science,  2014),  addressing  the  first  question.  We  are  now 
finalizing  a  manuscript  detailing  the  relationship  between  acinar 
and  ductal  adenocarcinoma.  We  find  that  they  are  generally 
clonally  related,  indicating  they  arise  from  the  same  cell  of  origin 
(Figure  1).  However,  there  are  specific  alterations  that  are  found 
much  more  frequently  in  ductal  foci  than  acinar  foci  and  likely 
contribute  to  the  histologic  appearance. 

While  addressing  this  aim  it  was  noted  that  mutations  identified  in 
metastases  were  found  at  low  frequency  in  high  grade,  localized 
foci  in  the  prostate.  In  collaboration  with  prostate  cancer  geneticists 
at  RIKEN  and  University  of  Tokyo  we  characterized  the  genomic 
heterogeneity  within  individual  index  foci  of  aggressive,  localized 
or  locally  advanced,  prostate  cancer.  Using  a  novel  tissue  banking 
method  (Gillard  et  al,  American  Journal  of  Translational  Research, 

2015)  we  performed  multiregion  genomic  analyisis  for  10  cases  of 
potentially  lethal  prostate  cancer.  Our  results  demonstrate  marked 
intratumoral  genetic  heterogeneity.  In  subclonal  populations  we 
identify  alterations  that  are  more  common  in  advanced,  castrate 
resistant  prostate  cancer.  Moreover,  there  is  heterogeneity  in 
genomic  biomarkers  which  predict  response  to  targeted  therapies 
(Figure  2).  The  manuscript  from  this  work  is  currently  under 
review. 


Figure  1.  Mutations  in 
coincident  ductal  and  acinar 
adenocarcinoma  of  the 
prostate.  Mutations  are 
represented  by  rows.  Darker 
green  represents  higher 
variant  allele  frequency. 
A=acinar,  D=ductal, 
M=metastasis. 
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Aim  2:  Characterize  the  molecular  relationship  between  metastatic  disease  and  primary  prostate  foci. 


This  aim  evaluated  the  relationship  between 
metastatic  and  primary  disease  in  two  scenarios. 
The  first  was  primary  prostate  cancer  and  nodal 
metastases,  where  we  found  a  close  relationship 
between  high  grade  disease  and  metastases 
(VanderWeele  et  al,  Cancer  Science,  2014),  and  in 
the  manuscript  currently  under  review. 


The  second  scenario  was  comparing  treatment- 
naive  disease,  late  stage  CRPC,  and  circulating 
tumor  cells  (CTCs).  Here  we  found  that  there  is 
progression  of  the  tumor  genome,  and  this  is 
captured  by  the  genomics  of  the  CTCs 
(Lack  et  al,  Journal  of  Translational  Medicine, 
2017). 


PROGRESS  REPORT  BY  SOW  TASKS 

Task  1:  Obtain  tissue  from  matched  prostate  cancer  foci  and  prepare  sequencing  libraries 

1 . 1  Completed,  manuscript  in  preparation 

1.2  Completed,  manuscript  published  (Lack  et  al,  JTM,  2017) 

1.3  Completed,  manuscript  in  preparation 

1.4  Completed,  manuscript  published  (Lack  et  al,  JTM,  2017) 

1 .5  Completed,  manuscript  under  review 

1 .6  In  progress 

Task  2.  Generate  next-generation  sequencing  data  on  Illumina  platform  with  massively  parallel  short 
(100  bp)  reads,  confirm  mutations  with  Sanger  sequencing  or  targeted  resequencing 

2.1  Completed,  manuscript  published  (VanderWeele  et  al,  Cancer  Science,  2014) 

2.2  Completed,  manuscript  published  (Lack  et  al,  JTM,  2017) 

2.3  Completed,  manuscript  published  (VanderWeele  et  al,  Cancer  Science,  2014) 

2.4  Completed,  manuscript  published  (Lack  et  al,  JTM,  2017)  and  manuscript  in  preparation 

2.5  Completed,  manuscript  under  review 

2.6  Completed,  manuscript  under  review 

2.7  Not  yet  completed 

Task  3.  Alignment  of  sequencing  reads,  calling  mutations,  and  constructing  phylogeny 

3.1  Completed,  manuscript  published  (VanderWeele  et  al,  Cancer  Science,  2014) 

3.2  Completed,  manuscript  published  (Lack  et  al,  JTM,  2017) 

3.3  Completed,  manuscript  published  (Lack  et  al,  JTM,  2017)  and  manuscript  in  preparation 

3.4  Completed,  manuscript  published  (Lack  et  al,  JTM,  2017)  and  manuscript  in  preparation 

3.5  Not  yet  completed 


Figure  2.  Heterogeneity  of  genomic  biomarkers 
within  the  index  focus  of  potentially  lethal, 
locally  advanced  prostate  cancer.  Cases  are  in 
rows,  biomarker  genes  in  columns. 
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Mentorship 


Dr.  Kathy  Kelly  is  my  primary  mentor,  and  she  has  continued  to  promote  my  professional 
development.  I  meet  with  her  individually  and  with  my  fellow  investigators  in  our  department  on  a 
monthly  basis.  The  most  recent  directions  of  my  research  have  been  strongly  positively  influenced  by 
her  work  and  the  resources  available  within  our  department. 


4  KEY  RESEARCH  ACCOMPLISHMENTS  FOR  YEAR  3: 

•  Acceptance  of  2  manuscripts,  including  one  on  evolution  of  the  genome  from  untreated 
disease  to  metastatic  disease  and  circulating  tumor  cells  (Aim  2b) 

•  Young  Investigator  Award  from  the  Prostate  Cancer  Foundation 

5.  CONCLUSION: 

Prostate  cancer  has  a  widely  varying  prognosis,  with  widely  varying  underlying  genomics.  The 
aggressive  variant  ductal  adenocarcinoma,  however,  does  not  appear  to  have  distinctive  molecular 
features  that  distinguish  it  from  conventional  acinar  adenocarcinoma,  including  similar  level  of 
Androgen  Receptor  activity.  The  similar  genomic  features  would  suggest  that  it  would  respond  to 
standard  therapies  in  a  manner  similar  to  aggressive  conventional  prostate  cancer. 

High  risk  prostate  cancer,  even  when  untreated,  harbors  alterations  similar  to  that  of  advanced, 
castrate  resistant  disease.  These  alterations  include  biomarkers  predictive  of  response  to  targeted 
therapy.  Testing  for  these  biomarkers  should  be  undertaken  with  caution,  however,  since  they  are 
often  subclonal.  This  implies  that  the  results  of  biomarker  tests  would  depend  on  where  the  tissue 
is  sampled,  and  the  indicated  therapy  may  not  be  effective  against  all  tumor  cells. 

Manuscripts  are  being  prepared  for  submission  to  report  on  the  findings  above.  Future  lab  work 
will  focus  on  in  vitro  organoid  models  and  in  vivo  xenograft  models  to  characterize  the  evolution 
of  tumor  cells  in  response  to  therapy. 


6.  PUBLICATIONS,  ABSTRACTS,  AND  PRESENTATIONS: 
a.  Peer-reviewed  Scientific  Journals 


Year  1 

1.  VanderWeele  DJ,  Brown  CB,  Taxy  JB,  Stadler  WM,  White  KP.  Low  grade  prostate  cancer  foci 
diverge  early  from  high  grade  and  metastatic  disease.  Cancer  Science.  2014  Aug;105(8):  1079-85. 

2.  Reyes  E*,  VanderWeele  DJ*.  Duggan  R,  Campanile  A,  Stadler  WM,  Vander  Griend  DJ, 
Szmulewitz  R.  Quantitative  characterization  of  androgen  receptor  expression  in  circulating 
tumor  cells  from  patients  with  metastatic  castrate  resistant  prostate  cancer.  Journal  of 
Translational  Medicine.  2014  Nov  26;12(1):313. 

Year  2 
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1.  Packiam  V,  Patel  SG,  Pariser  JJ,  Richards  KA,  Weiner  AB,  Paner  GP,  VanderWeele  DJ,  Zagaja 
GP,  Eggener  SE.  “Contemporary  Population-Based  Comparison  of  Localized  Ductal 
Adenocarcinoma  and  High-Risk  Acinar  Adenocarcinoma  of  the  Prostate”.  Urology.  2015  Jul 
18. 

2.  VanderWeele  DJ.  Paner  GP,  Fleming  GF,  Szmulewitz  RZ.  Sustained  complete  response  to 
cytotoxic  therapy  and  the  PARP  inhibitor  veliparib  in  metastatic  castration  resistant  prostate 
cancer  -  a  case  report.  Frontiers  in  Oncology.  2015  5:169. 

3.  Gillard  M,  Tom  W,  Antic  T,  Paner  GP,  Lingen  M,  VanderWeele  DJ.  “Next-gen  tissue: 
preservation  of  molecular  and  morphological  fidelity  in  prostatectomy  specimens”.  American 
Journal  of  Translational  Research.  2015  Jul  15;  7(7):  1227-1235. 

4.  McCann  S,  Jiang  Y,  Fan  X,  Wang  J,  Antic  T,  VanderWeele  DJ,  Oto  A.  “Quantitative  Multi- 
parametric  MRI  Features  and  PTEN  Expression  in  Prostate  Cancer”.  American  Journal  of 
Roentgenology.  2016  Mar;206(3):559-65. 

Year  3 

1.  LackJ,  Gillard  M,  CamM,  Paner  GP,  VanderWeele-DJ.  “Circulating  tumor  cells  capture  disease 
evolution  in  advanced  prostate  cancer”.  Journal  of  Translational  Medicine.  2017  Feb  23;15(1)44 


b.  Reviews  and  Book  Chapters 


Year  2 

1.  VanderWeele  DJ,  Stadler  WM.  “The  time  for  molecular  phenotyping  of  metastatic  castrate- 
resistant  prostate  cancer  is  now”.  Clinical  Advances  in  Hematology  &  Oncology.  2015  13(9),  p. 
580. 

Year  3 

1 .  VanderWeele  DJ,  Turkbey  B,  Sowalsky  A.  “Precision  management  of  localized  prostate  cancer”. 
Expert  Review  of  Precision  Medicine  and  Drug  Development.  2016  Dec;l(6):505-15. 

2.  VanderWeele  DJ,  McNerney  M,  LeBeau  M.  “Chromosomal  Aberrations  in  Cancer”.  Holland- 
Frei  Cancer  Medicine,  9th  edition.  Bast  RC,  Croce  CM,  Hait  WN,  et  al,  editors. 


c.  Presentations 


Year  3 

1.  VanderWeele  DJ,  Finney  R,  Katayama  K,  Gillard  M,  Paner  G,  Imoto  S,  Yamaguchi  R,  Wheeler 
D,  Cam  M,  Maejima  K,  Sasaki-Oku  A,  Nakano  K,  Tanaka  H,  Pontier  A,  Grigoryev  D,  Kubo  M, 
Ratain  M,  Miyano  S,  Nakagawa  H.  “Local  Enrichmnet  of  Multiple  Clinically  Significant 
Genetic  Alterations  within  Index  Foci  of  Locally  Advanced  Porstate  Cancer.”  DOD  IMPaCT 
meeting.  Aug.,  2016.  Oral  presentation. 

2.  VanderWeele  DJ,  Finney  R,  Katayama  K,  Gillard  M,  Paner  G,  Imoto  S,  Yamaguchi  R,  Wheeler 
D,  Cam  M,  Maejima  K,  Sasaki-Oku  A,  Nakano  K,  Tanaka  H,  Pontier  A,  Grigoryev  D,  Kubo  M, 
Ratain  M,  Miyano  S,  Nakagawa  H.  “Treatment-naive,  locally  advanced  prostate  cancer  contains 
numerous  subclonal  clinically  significant  alterations.”  Prostate  Cancer  Foundation  annual 
scientific  retreat.  Oct.,  2016.  Poster  presentation. 
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3.  Vander  Weele  DJ,  Lack  J,  Gillard  M,  Cam  M,  Paner  GP.  “Genomic  analysis  of  untreated 

prostate  adenocarcinoma,  treatment-emergent  neuroendocrine  carcinoma,  and  circulating  tumor 
cells  (CTCs).”  Genitourinary  Cancer  Symposium.  Feb.,  2017.  Poster  presentation. 


7.  INVENTIONS,  PATENTS  AND  LICENSES: 

Nothing  to  report 

8.  REPORTABLE  OUTCOMES: 

Nothing  to  report 

9.  OTHER  ACHIEVEMENTS: 

The  data  from  this  work,  and  the  protected  time  insured  from  this  award,  has  led  to  the  successful 
application  for  a  Young  Investigator  Award  from  the  Prostate  Cancer  Foundation. 

10.  OPPORTUNITIES  FOR  TRAINING: 

As  per  the  SOW,  I  have  completed  the  following  courses 

•  Responsible  Conduct  of  Research 

•  ECEV  35600  Principles  of  Population  Genetics  I 

•  ECEV  35901  Evolutionary  Genomics 

•  Fundamentals  of  Clinical  Research 

•  HGEN  47400  Introduction  to  Probability  and  Statistics  for  Geneticists 

•  STAT  22000  Statistical  Methods  &  Their  Applications 

•  Introduction  to  the  Principles  and  Practice  of  Clinical  Research 

•  Ethical  and  Regulatory  Aspects  of  Clinical  Research 

I  have  yet  to  complete 

•  Principles  of  Clinical  Pharmacology 


11.  OPPORTUNITIES  FOR  PROFESSIONAL  DEVELOPMENT: 

I  had  the  opportunity  to  attend  and  participate  in  international  conferences  during  the  third  year  of  this 
award.  I  gave  an  oral  presentation  at  the  Department  of  Defense  IMPaCT  meeting  in  August,  2016.  I 
received  valuable  feedback  on  the  project  after  this  presentation. 

I  gave  poster  presentations  at  the  Prostate  Cancer  Foundation  annual  scientific  retreat  and  at  the  GU 
Cancers  Symposium  (GU  ASCO).  The  work  presented  at  GU  ASCO  was  published  shortly 
thereafter,  describing  the  relationship  between  treatment-naive  prostate  cancer,  metastatic  disease,  and 
CTCs.  Both  of  these  meetings  have  clinical  and  translational  components. 

Though  not  directly  related  to  my  work  in  the  lab,  I  gave  an  oral  presentation  on  prostate  cancer 
genetics  and  its  clinical  implications  at  the  International  Prostate  Imaging  Symposium.  The  work 
supported  by  this  award  led  to  the  recognition  of  my  expertise  in  this  area  and  the  invitation  to  speak. 
Similarly,  the  expertise  gained  from  this  award  supported  my  successful  application  for  a  Young 
Investigator  Award  from  the  Prostate  Cancer  Foundation. 
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Finally,  I  have  established  collaborations  across  diverse  fields  to  conduct  a  pre-clinical  trial  evaluating 
a  novel  class  of  Androgen  Receptor  (AR)  antagonist,  and  to  conduct  correlative  studies  in  conjunction 
with  a  clinical  trial  of  neoadjuvant  AR-targeted  therapy. 


12.  REFERENCES: 

1 .  VanderWeele  DJ,  Brown  CB,  Taxy  JB,  Stadler  WM,  White  KP.  Low  grade  prostate  cancer  foci 
diverge  early  from  high  grade  and  metastatic  disease.  Cancer  Science.  2014  Aug;105(8):  1079-85. 

2.  Gillard  M,  Tom  W,  Antic  T,  Paner  GP,  Lingen  M,  VanderWeele  DJ.  “Next-gen  tissue: 
preservation  of  molecular  and  morphological  fidelity  in  prostatectomy  specimens”.  American 
Journal  of  Translational  Research.  2015  Jul  15;  7(7):1227-1235. 

3.  LackJ,  Gillard  M,  CamM,  Paner  GP,  VanderWeele-DJ.  “Circulating  tumor  cells  capture  disease 
evolution  in  advanced  prostate  cancer”.  Journal  of  Translational  Medicine.  2017  Feb  23;15(1)44. 


13  APPENDICES: 


1 .  VanderWeele  DJ,  Brown  CB,  Taxy  JB,  Stadler  WM,  White  KP.  Low  grade  prostate  cancer  foci 
diverge  early  from  high  grade  and  metastatic  disease.  Cancer  Science.  2014  Aug;105(8):  1079-85. 

2.  Gillard  M,  Tom  W,  Antic  T,  Paner  GP,  Lingen  M,  VanderWeele  DJ.  “Next-gen  tissue: 
preservation  of  molecular  and  morphological  fidelity  in  prostatectomy  specimens”.  American 
Journal  of  Translational  Research.  2015  Jul  15;  7(7):  1227-1235. 

3.  LackJ,  Gillard  M,  CamM,  Paner  GP,  VanderWeele-DJ.  “Circulating  tumor  cells  capture  disease 
evolution  in  advanced  prostate  cancer”.  Journal  of  Translational  Medicine.  2017  Feb  23;15(1)44. 

4.  VanderWeele  DJ  et  al,  manuscript  under  review  at  JNCI. 
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Understanding  the  developmental  relationship  between  indolent  and  aggressive 
tumors  is  central  to  understanding  disease  progression  and  making  treatment 
decisions.  For  example,  most  men  diagnosed  with  prostate  cancer  have  clinically 
indolent  disease  and  die  from  other  causes.  Overtreatment  of  prostate  cancer 
remains  a  concern.  Here  we  use  laser  microdissection  followed  by  exome 
sequencing  of  low-  and  high-grade  prostate  cancer  foci  from  four  subjects,  and 
metastatic  disease  from  two  of  those  subjects,  to  evaluate  the  molecular  relation¬ 
ship  of  coincident  cancer  foci.  Seventy  of  79  (87%)  high-confidence  somatic  muta¬ 
tions  in  low-grade  disease  were  private  to  low-grade  foci.  In  contrast,  high-grade 
foci  and  metastases  harbored  many  of  the  same  mutations.  In  cases  in  which 
there  was  a  metastatic  focus,  15  of  80  (19%)  high-confidence  somatic  mutations 
in  high-grade  foci  were  private.  Seven  of  the  80  (9%)  were  shared  with  low- 
grade  foci  and  65  (82%)  were  shared  with  metastatic  foci.  Notably,  mutations  in 
cancer-associated  genes  and  the  p53  signaling  pathway  were  found  exclusively 
in  high-grade  foci  and  metastases.  The  pattern  of  mutations  is  consistent  with 
early  divergence  between  low-  and  high-grade  foci  and  late  divergence  between 
high-grade  foci  and  metastases.  These  data  provide  insights  into  the  develop¬ 
ment  of  high-grade  and  metastatic  prostate  cancer. 
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Cancers  are  thought  to  progress  from  low  grade  to  higher 
grade/ leading  to  efforts  to  detect  and  eradicate  low- 
grade  disease  in  order  to  prevent  morbidity  and  mortality  from 
high-grade  disease.  The  grade  of  prostate  cancer  is  character¬ 
ized  by  the  Gleason  score  on  a  scale  of  2-10,  with  10  repre¬ 
senting  the  most  poorly  differentiated  tumors.  Although 
initially  described  as  intermediate  grade,  Gleason  6  cancers  are 
clinically  indolent  and  have  an  overwhelmingly  favorable 
prognosis. (2,3)  For  the  purposes  of  this  report,  Gleason  6  foci 
are  regarded  as  low  grade. 

Approximately  238  590  men  were  diagnosed  with  prostate 
cancer  in  the  US  in  2013,  the  majority  of  them  with  low-grade 
disease.  It  is  estimated  that  29  720  men  will  die  from  prostate 
cancer,  which  is  among  the  lowest  case/fatality  ratio  of  any 
cancer. (4)  Although  prostate  specific  antigen  (PSA)  screening  is 
able  to  identify  early  stage  prostate  cancer,  most  men  diag¬ 
nosed  with  prostate  cancer  die  from  competing  causes.  It  is 
therefore  unclear  whether  early  detection  associated  with  PSA 
screening  decreases  mortality. (5,6)  Current  clinical  treatment 
guidelines  support  active  surveillance  for  patients  with 
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low-grade  prostate  cancer  for  which  the  risk  of  tumor  progres¬ 
sion  is  low.(7)  Recently  some  guidelines  have  even  recom¬ 
mended  against  PSA  screening. (8)  Failure  to  follow  these 
guidelines  is  likely  due,  in  part,  to  concern  that  low-grade  foci 
progress  to  higher  grade  and  metastatic  disease. 

Several  studies  have  contributed  to  a  growing  catalog  of 
genetic  variation  in  prostate  cancer, (9_14)  giving  insight  into 
single  nucleotide  variations,  rearrangements  and  copy  number 
variation  prevalent  in  this  disease.  Aside  from  common  rear¬ 
rangements  involving  the  ETS  family  of  transcription  factors, 
especially  v-ets  avian  erythroblastosis  virus  E26  oncogene 
homolog  (ERG),  the  mutational  spectrum  of  prostate  cancers  is 
diverse. (10)  Although  increasing  attention  is  being  paid  to  the 
overall  molecular  heterogeneity  of  cancers  within  a  single 
patient,  there  is  limited  data  in  this  regard  specifically  for  pros¬ 
tate  cancer.  This  is  especially  relevant  for  prostate  cancer, 
which  is  often  multifocal. 

More  specifically,  it  is  not  clear  if  prostate  cancer  develops 
as  a  multifocal  disease  in  which  the  foci  are  clonally  related 
and  share  the  vast  majority  of  somatic  mutations,  or  if  it  is  a 
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multicentric  disease,  with  each  focus  having  its  own  set  of  pri¬ 
vate  somatic  mutations  and  each  representing  an  independent 
primary  lesion.  Recent  studies  have  examined  multifocal 
disease  and  found  evidence  of  a  monoclonal  origin  of  coinci¬ 
dent  foci,(15)  even  when  the  foci  have  different  Gleason  pat¬ 
terns. (16,17)  This  conclusion  is  somewhat  surprising,  given  the 
difference  in  clinical  behavior  of  different  Gleason  patterns. 

The  present  study  targets  multifocal  disease,  specifically 
with  low-grade  (Gleason  6)  and  high-grade  (Gleason  8  or 
higher)  foci,  and  compares  the  molecular  relationship  of  these 
foci  with  each  other  and  with  synchronous  metastatic  disease. 
To  expand  the  number  of  specimens  available  for  the  present 
study,  prostatectomy  samples  that  had  been  formalin  fixed  and 
paraffin  embedded  (FFPE)  were  used,  with  modifications  made 
to  standard  protocols  to  increase  the  yield  of  material  used  for 
sequencing.  Exome  sequencing  was  used  to  achieve  a  compre¬ 
hensive  evaluation  of  all  genes  harboring  somatic  point  muta¬ 
tions.  The  findings  indicate  that  there  is  early  divergence  or 
complete  independence  of  low-grade  and  high-grade  disease, 
but  late  divergence  of  high-grade  disease  and  synchronous 
metastases. 

Materials  and  Methods 

Tissue  and  laser  capture  microdissection  (LMD).  Paraffin- 
embedded  prostate  and  lymph  node  tissue  that  was  1-5  years 
old  and  matched  FFPE  and  frozen  normal  prostate  tissue  was 
obtained  from  the  Human  Tissue  Resource  Center  (HTRC)  at 
the  University  of  Chicago  with  the  approval  of  the  Internal 
Review  Board.  Tissue  was  reviewed  by  a  genitourinary  pathol¬ 
ogist  (J.B.T.)  and  regions  of  interest  were  marked  for  further 
evaluation.  Four  specimens  had  sufficient  low-  and  high-grade 
disease  for  LMD.  Two  specimens  (PrCa  6  and  PrCa  18)  also 
had  sufficient  lymph  node  metastatic  disease.  Tumor  foci  and 
histologically  normal  prostate  were  microdissected  from  10- 
pm  sections  using  a  Leica  Laser  Microdissection  (Wetzlar, 
Germany)  system. 

Exome  library  preparation.  Genomic  FFPE  DNA  was  isolated 
using  a  QIAamp  DNA  FFPE  Tissue  Kit  (Qiagen,  Valencia, 
CA,  USA)  following  the  manufacturer’s  instructions,  but 
increasing  the  incubation  in  Proteinase  K  to  overnight  and 
reducing  the  90°C  incubation  to  30  min.  DNA  from  frozen  tis¬ 
sue  was  isolated  using  a  QIAamp  Mini  Kit  (Qiagen).  DNA 


Table  1.  Subject  characteristics 


Case 

Patient  age  (years) 

Race 

Stage 

Gleason  -  glandf 

Library 

Gleason  -  focusj 

PrCa  18 

65 

Caucasian 

pT3bN1  Mx 

4  +  3 

Normal 

NA 

Low  grade 

6 

High  grade 

8 

Metastasis 

NA 

PrCa  14 

54 

African- 

pT3bN1  Mx 

4  +  3 

Normal 

NA 

American 

Low  grade 

6 

High  grade 

8 

PrCa  6 

70 

Caucasian 

pT3bN1  Mx 

4  +  3  (5) 

Normal 

NA 

Low  grade 

6 

High  grade 

8 

Metastasis 

NA 

PrCa  25 

59 

Asian 

pT2cN0Mx 

3  +  4 

Normal 

NA 

Low  grade  6 

High  grade  9 


tGleason  score  assigned  on  pathological  review  of  the  prostatectomy  specimen.  JGleason  score  of  the  focus  isolated  by  microdissection.  PrCa  6 
had  a  Gleason  score  of  4  +  3  with  a  minor  component  of  5.  PrCa  14  had  metastatic  disease  in  a  lymph  node,  which  was  not  of  sufficient  mass  to 
be  isolated  for  sequencing.  NA,  not  applicable. 


was  fragmented  using  a  Covaris  S2  sonicator  (Woburn,  MA, 
USA),  with  settings  of  10%  duty  cycle,  intensity  4,  200  cycles 
per  burst  and  120  s.  The  fragmented  DNA  was  blunt-end 
repaired  using  T4  Polymerase,  Klenow  and  T4  PNK,  and  ade- 
nylated.  Fragments  were  ligated  to  adaptors  ISl_adapter_P5.F 
and  IS3_adapter_P5  +  P7.R(18)  and  an  index  tag  added.  The 
product  was  amplified  and  cleaned  up  using  Agencourt  AM- 
Pure  XP  beads  (Beckman  Coulter,  Indianapolis,  IN,  USA). 

Libraries  were  run  on  a  BioAnalyzer  (Agilent,  Santa  Clara, 
CA,  USA)  to  confirm  fragment  size.  The  libraries  from  the 
metastases  did  not  undergo  further  size  selection.  The  remain¬ 
ing  libraries  underwent  selection  using  an  E-gel  (Invitrogen, 
Grand  Island,  NY,  USA)  followed  by  PCR  amplification. 

Enrichment  for  exomes  was  performed  using  NimbleGen 
SeqCap  EZExome  v2.0  (Roche,  Madison,  WI,  USA)  following 
the  manufacturer’s  instructions.  The  exome-enriched  library 
DNA  was  amplified  and  run  on  a  BioAnalyzer  to  confirm  frag¬ 
ment  size.  The  resulting  libraries  were  run  on  a  HiSeq  2000 
(Illumina,  San  Diego,  CA,  USA)  in  a  highly  multiplexed 
1  x  50  run  and  a  subsequent  2  x  100  run. 

Sequencing  analysis  and  identification  of  somatic  mutations. 
Overlapping  or  adapter-containing  reads  were  clipped  using 
SeqPrep-b83fd00.  Reads  were  quality  trimmed  and  aligned  to 
hgl8  using  bwa  0.5.9(19)  and  merged  with  duplicates  removed 
using  Samtools  0.1.18.(2O)  Genome  Analysis  Toolkit(21)  was 
used  for  local  realignment  and  quality  score  recalibration. 
Somatic  mutations  were  identified  using  VarScan  v2.2.8.(22) 
High-confidence  somatic  mutations  were  defined  as  variants 
with  coverage  of  at  least  10  x,  no  supporting  reads  in  normal, 
at  least  four  supporting  reads  in  the  tumor  and  with  supporting 
reads  on  both  strands.  Germline  single-nucleotide  polymor¬ 
phisms  (SNPs)  identified  in  dbSNP135,  1000  genomes  and  ex¬ 
omes  of  normal  controls  published  by  the  Max  Planck  Institute 
(MPI)(23-25)  were  removed.  SAMtools  was  used  to  determine 
the  number  of  total  and  variant  reads  in  matched  libraries  for 
positions  mutated  in  low-grade  foci.  Sample  ethnicity  was 
determined  by  principle  component  analysis  using  common, 
germline,  coding  variants  shared  with  samples  from  the  1000 
genomes  project. (24) 

Custom  capture  and  identification  of  previously  identified 
mutations.  DNA  from  seven  additional  foci  from  PrCa  6 
were  isolated  using  LMD  (Leica  Microsystems  LMD6500). 
Custom  capture  probes  targeting  high-confidence  mutations 
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identified  through  exome  sequencing  were  obtained  from 
Roche.  Sequencing  libraries  from  the  seven  loci  were  pre¬ 
pared  using  the  High-Throughput  Genome  Analysis  Core  at 
the  University  of  Chicago,  pooled,  captured  following  the 
manufacturer’s  instructions,  and  sequenced  on  one  lane  of  a 
HiSeq2000  (Illumina).  Reads  were  analyzed  with  SeqPrep- 
b83fdOO,  aligned  to  hgl8  using  bwa,  and  merged  with 
duplicates  removed  using  Samtools  0.1.18.  Samtools  was 
used  to  determine  the  number  of  total  and  variant  reads  in 
each  library  for  each  variant  position.  Variants  were  consid¬ 
ered  present  if  >10%  of  reads  supported  the  variant. 

Confirmation  of  somatic  mutations  and  immunohistochemisty 
(IHC).  Mutations  were  visualized  using  IGV  2.0.(26)  Genomic 
DNA  from  each  focus  underwent  amplification  using  Genome  - 
Plex  WGA2  (Sigma,  St.  Louis,  MO,  USA)  and  primers  target¬ 
ing  a  subset  of  variants  were  designed  using  Primer3Plus  and 
obtained  from  IDT  (Coral  ville,  I  A,  USA).  The  PCR  products 
were  sequenced  in  the  University  of  Chicago  DNA  Sequencing 
and  Genotyping  Facility.  ERG  IHC  was  performed  at  the 


HTRC  using  anti-ERG  antibody  EPR3864  (Abeam, 
Cambridge,  MA,  USA). 

Results 

To  evaluate  the  molecular  relationship  between  low-  and  high- 
grade  prostate  cancer,  we  identified  multifocal  formalin-fixed, 
paraffin-embedded  specimens  with  both  low-grade  and  high- 
grade  disease  with  sufficient  mass  for  laser  capture  microdis¬ 
section.  The  subjects  from  whom  the  specimens  came  were 
representative  of  patients  with  high-risk  localized  prostate  can¬ 
cer  (Table  1). 

Exome  sequencing  was  performed  on  matched  cancer  foci 
and  histologically  normal  prostate  glands  from  four  specimens 
(Fig.  1).  To  improve  the  yield  of  DNA  for  downstream 
sequencing,  proteinase  K  digestion  was  increased  to  18  h,  90° 
incubation  was  reduced  to  30  min  and  size  selection  by  gel 
extraction  was  replaced  with  size  selection  with  an  E-gel  or 
eliminated  altogether.  On  average,  88%  of  RefSeq  coding 


Fig.  1.  Histology  of  coincident  prostate  cancer  foci.  Representative  H&E  stains  and  illustrations  representing  the  prostate  cancer  foci  that  were 
laser  microdissected.  Two  cancer  foci  and  uninvolved  prostate  glands  were  isolated  from  PrCa  18  (a),  PrCa  14  (b),  PrCa  6  (c)  and  PrCa  25  (d).  In 
addition,  a  metastatic  (Met)  focus  was  isolated  from  PrCa  18  and  PrCa  6.  For  PrCa  14  (b),  the  foci  were  from  different  levels  of  the  prostate.  For 
the  other  three  specimens  the  foci  were  at  the  same  level.  Light  gray,  histologically  normal  prostate;  dark  gray,  low-grade  cancer  focus;  striped, 
high-grade  cancer  focus;  checked,  metastatic  focus  from  a  lymph  node  removed  at  the  time  of  prostatectomy. 
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Fig.  2.  Mutation  characteristics,  (a)  Number  of 
high-confidence  somatic  mutations  across  all  foci. 
Non-silent,  non-silent  mutations;  Unique,  number 
of  unique  genes  harboring  a  non-silent  mutation; 
Reported,  gene  reported  to  be  mutated  in 
references  9-12  and  14.  (b)  Spectrum  of  unique 
high  confidence  somatic  mutations  across  all  foci, 
(c)  Number  of  high-confidence  somatic  mutation 
types  in  all  prostate  samples.  L,  low-grade  focus;  H, 
high-grade  focus;  M,  metastatic  focus;  NS, 
nonsynonymous;  SYN,  synonymous;  UTR3,  3' 
untranslated  region;  UTR5,  5'  untranslated  region; 
SJ,  splice  junction;  STOP,  premature  stop  codon,  (d) 
Confirmation  of  high-confidence  somatic  mutations 
with  capillary  sequencing  by  number  of  total  and 
variant  reads  (gray,  not  confirmed;  black, 
confirmed).  Confirmation  rate,  0.85. 


bases  were  covered  at  10  x  or  greater.  A  mean  of  30.1  somatic 
coding  mutations  was  identified  per  cancer  focus  (Supporting 
Information  Table  SI). 

To  evaluate  potential  bias  introduced  by  formalin  fixation, 
matched  FFPE  and  frozen  (Optimal  Cutting  Temperature 
compound  [OCT] -embedded)  samples  of  normal  prostate  were 
subjected  to  exome  sequencing.  The  number  of  false  positive 
variant  reads  was  similar  (Fig.  SI  A),  as  were  the  types  of 
mutations  (Fig.  SIB),  suggesting  that  there  was  little  bias 
introduced  from  formalin  fixation. 

Three  hundred  and  one  somatic  mutations  were  identified 
across  10  foci  (Table  S2).  Of  these,  71  were  synonymous  and 
unlikely  to  alter  gene  function  (Fig.  2a).  Of  174  unique  genes 
harboring  somatic  mutations  that  were  likely  to  be  functional,  73 
were  not  reported  in  five  previous  reports  of  mutations 
in  localized  and  advanced  disease/9-12,1  }  The  majority  of 
variants  identified  were  transitions  leading  to  nonsynonymous 
mutations  (Fig.  2b, c),  as  expected.  A  subset  of  20  high- 
confidence  somatic  mutations  was  confirmed  in  tumor  and  unin¬ 
volved  tissue  with  capillary  sequencing  (Fig.  2d).  Those  that 
were  not  confirmed  by  capillary  sequencing  had  low  variant 
allele  frequency  and  were  approximately  the  level  of  detection 
for  capillary  sequencing  (-20%).  Nevertheless,  the  mutation 
(0.82/Mb)  and  validation  rate  (85%)  using  fixed  prostate 
specimens  in  the  present  study  are  similar  to  previous  reports 
(0.9/Mb  and  91%,  respectively)  using  frozen  tissue. (9)  The  high 
fidelity  of  this  sequencing  supports  the  use  of  formalin-fixed  tis¬ 
sue  for  next-generation  sequencing  studies. 

The  vast  majority  of  somatic  mutations  in  low-  and  high-grade 
foci  were  private  to  the  focus  in  which  they  were  identified.  Nine 
of  79  (1 1%)  high-confidence  somatic  mutations  in  the  low-grade 
foci  were  shared  with  a  high-grade  focus,  and  nine  of  139  (7%) 
high-confidence  somatic  mutations  in  the  high-grade  foci  were 
shared  with  a  low-grade  focus  (Fig.  3).  Capillary  sequencing 
confirmed  that  these  nine  shared  somatic  mutations  were  not 
germline  variants.  In  two  subjects  (14  and  18)  all  of  the  muta¬ 
tions  were  private.  In  the  other  two  subjects  (6  and  25)  a  small 


minority  of  mutations  was  shared.  In  contrast  to  the  low-grade 
foci,  67  of  80  (84%)  high-confidence  somatic  mutations  identi¬ 
fied  in  high-grade  foci  were  shared  with  metastatic  disease  and 
metastatic  foci  had  few  mutations  beyond  those  identified  in 
high-grade  disease  (Fig.  3a,c). 

To  evaluate  more  broadly  the  clonal  relationship  of  multifo¬ 
cal  disease,  custom  capture  baits  targeting  variants  identified 
through  exome  sequencing  of  PrCa  6  were  used  to  capture  the 
DNA  from  seven  additional  cancer  foci  from  that  subject 
(Table  S3).  The  average  coverage  depth  of  these  variants  was 
208  and  78  previously  identified  variants  were  covered  at 
>70  x  in  each  focus.  Three  high-grade  foci  harbored  a  median 
of  56  previously  identified  variants,  three  intermediate-grade 
foci  harbored  a  median  of  15  variants  and  one  low-grade  focus 
harbored  three  variants.  As  in  the  whole  exome  data  a  small 
minority  of  mutations  was  shared,  with  just  two  of  78  (2.6%) 
variants  found  to  be  ubiquitous  among  the  seven  cancer  foci. 
Fourteen  of  78  (17.9%)  were  not  identified  in  any  of  the  seven 
additional  foci  (Fig.  S2).  Eight  of  these  14  were  among  the  12 
variants  originally  identified  as  private  to  the  low-grade  focus. 

Rearrangements  involving  ERG  are  the  most  common  genetic 
derangement  in  prostate  cancer  and  are  thought  to  occur  early  in 
tumorigenesis.  These  events  are  not  typically  captured  by  exome 
sequencing.  Therefore  we  performed  IHC  for  ERG  to  look  for 
overexpression,  which  is  typically  restricted  to  cells  harboring 
an  ERG  rearrangement.  Of  the  two  specimens  exhibiting  posi¬ 
tive  ERG  staining,  PrCa  6  showed  consistently  positive  staining 
in  all  tumor  foci  (Fig.  4).  In  contrast,  in  PrCa  18  the  low-grade 
focus  was  negative,  the  high-grade  focus  was  equivocal  with 
few  cells  positive  for  ERG  and  the  metastatic  focus  was  more 
widely  focally  positive. 

Together  with  the  immunohistochemical  analysis,  the  pat¬ 
terns  of  shared  and  private  mutations  indicate  early  divergence 
of  low-  and  high-grade  foci.  In  two  cases  no  shared  mutations 
were  identified  and  independent  origin  of  low-  and  high-grade 
disease  cannot  be  ruled  out.  In  contrast,  there  is  late  diver¬ 
gence  of  high-grade  foci  and  metastatic  disease.  Consistent 
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Fig.  3.  The  molecular  relationship  of  coincident  foci.  Venn  diagrams 
(left)  depicting  the  pattern  of  shared  and  private  high-confidence 
somatic  mutations  and  phylogenetic  trees  (right)  depicting  the  rela¬ 
tionship  of  coincident  foci  for  PrCa  18  (a),  PrCa  14  (b),  PrCa  6  (c)  and 
PrCa  25  (d).  The  number  of  high-confidence  somatic  mutations  is 
labeled  within  each  Venn  diagram.  In  (c),  there  are  no  mutations  that 
are  shared  between  the  low-grade  focus  and  the  metastatic  focus  and 
also  not  shared  with  the  high-grade  focus.  Within  each  phylogenetic 
tree,  branch  length  is  proportional  to  the  number  of  mutations.  Gray, 
uninvolved  prostate;  blue,  low-grade  focus;  green,  high-grade  focus; 
purple,  metastatic  focus;  black,  theoretical  common  progenitor. 


with  this,  there  were  no  somatic  mutations  that  were  shared 
between  low  grade  and  metastatic  disease  and  also  not  shared 
with  high-grade  disease. 

Prostate  tumor  foci  consist  of  an  admixture  of  tumor  cells  and 
stromal  cells.  Despite  microdissection  of  the  tumor  foci,  the  var¬ 


iant  allele  frequencies  in  the  low-grade  foci  were  modestly  lower 
than  the  high-grade  lesions  (Fig.  S3),  suggesting  higher  contami¬ 
nation  with  non-tumor  cells.  The  contaminating  stromal  cells 
may  interfere  with  the  ability  to  call  somatic  mutations,  espe¬ 
cially  in  the  low-grade  foci,  leading  to  a  falsely  elevated  number 
of  apparently  private  mutations  and  obscuring  the  clonal  rela¬ 
tionship  among  concurrent  foci.  Therefore  for  each  position 
identified  as  a  private  high-confidence  somatic  mutation  in  the 
low-grade  foci,  the  number  of  sequencing  reads  supporting  that 
mutation  in  the  matched  high  grade  and  metastatic  foci  was 
determined,  relaxing  filters  based  on  DNA  strand  and  quality 
scores  (Fig.  S4).  Of  the  23  high-confidence  somatic  mutations  in 
PrCa  18  and  PrCa  6  identified  as  private  mutations,  none  had 
two  or  more  variant  reads  in  high  grade  or  metastatic  foci.  Thus, 
even  with  a  less  conservative  definition  of  mutation  (two  variant 
reads),  high-grade  foci  and  metastases  harbored  few  variants 
seen  in  low-grade  foci  and  the  majority  of  variants  in  the  low- 
grade  foci  was  private.  This  was  corroborated  with  evidence 
from  capillary  sequencing.  Three  mutations  private  to  low-grade 
foci  and  two  shared  with  other  foci  were  among  the  subset  of 
mutations  confirmed  by  capillary  sequencing. 

Two  genes  exhibited  recurrent  nonsynonymous  mutations  or 
mutations  in  the  untranslated  region  (Table  S4).  Low  density 
lipoprotein  (LDL)  receptor-related  protein  IB,  which  encodes 
an  LDL-family  receptor,  is  frequently  deleted  in  multiple 
malignancies. (27)  However,  recent  work  suggests  its  high  fre¬ 
quency  of  mutation  may  not  imply  a  significant  role  in  can¬ 
cer/28^  ZNF717  is  one  of  a  number  of  zinc  finger  proteins  and 
its  role  in  malignancy  is  unclear. 

Comparison  of  genes  with  nonsynonymous  mutations  in  our 
data  to  the  Catalogue  Of  Somatic  Mutations  In  Cancer  (COS¬ 
MIC)  cancer  gene  census  revealed  three  genes  in  common 
(Table  2).  These  genes  include  TP53,  the  most  frequently 
mutated  gene  in  advanced  cancers  including  prostate  cancer/1  } 
as  well  as  ATM  and  SS18L1.  All  of  these  three  genes  were 
mutated  in  high  grade  or  metastatic  disease,  but  not  in  low-grade 
disease.  Gene  ontology  analysis(29)  revealed  enrichment  for 
genes  in  the  p53  signaling  pathway  among  high-grade  foci 
(adjusted  P-value  0.037),  including  TP53  and  ATM.  Cell  cycle 
genes  were  also  enriched  (adjusted  P-value  0.054).  No  other 
pathways  were  statistically  significantly  enriched. 


Discussion 

Given  the  disparate  outcomes  between  men  with  low-grade  pros¬ 
tate  cancer  and  those  with  high-grade  disease,  a  fundamental 
question  in  prostate  cancer  management  is  the  possible  molecu¬ 
lar  relationship  among  grades  and  metastatic  disease.  Multiple 
lines  of  evidence  suggest  multifocal  disease  is  often  indepen¬ 
dent/30,3^  However,  recent  reports  have  demonstrated  shared 
molecular  derangements  of  coincident  low-  and  high-grade 
foci.(16,17)  The  present  study  identifies  mutations  shared  between 
low-  and  high-grade  cancer  foci,  supporting  the  ability  of  a  sin¬ 
gle  progenitor  to  give  rise  to  both  low-  and  high-grade  disease. 
However,  the  relationship  between  these  foci  was  distant  relative 
to  the  relationship  between  high-grade  disease  and  synchronous 
metastases.  In  fact,  in  two  of  the  four  cases  it  cannot  be  ruled  out 
that  the  foci  arose  from  independent  origins.  The  early  diver¬ 
gence  of  low-  and  high-grade  disease  may  help  explain  distinct 
differences  in  their  clinical  behavior.  Indeed,  we  found  no  evi¬ 
dence  of  direct  progression  from  low  grade  to  metastatic  disease. 
Instead,  there  was  an  overwhelming  fraction  of  shared  mutations 
between  high  grade  and  metastatic  disease,  which  is  consistent 
with  late  divergence  of  these  high-grade  foci  (Fig.  5). 
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Low  grade  High  grade  Metastasis 


Fig.  4.  ERG  expression.  Immunohistochemical 
analysis  for  ERG  in  coincident  foci  from  PrCa  6  and 
PrCa  18. 


Table  2.  Overlap  of  mutated  genes  with  the  COSMIC  cancer  gene  census 


Subject 

Histology 

Chromosome 

Hg18  position 

Reference 

Variant 

Symbol 

GERP 

Type 

aa  Change 

6 

High 

11 

107675775 

G 

A 

ATM 

5.22 

STOP 

W1710* 

6 

LN 

11 

107675775 

G 

A 

ATM 

5.22 

STOP 

W1710* 

6 

High 

20 

60169915 

C 

G 

SS18L1 

3.98 

NS 

P87R 

6 

LN 

20 

60169915 

C 

G 

SS18L1 

3.98 

NS 

P87R 

14 

High 

17 

7517833 

C 

A 

TP53 

4.78 

NS 

C145F 

Three  genes  harboring  a  nonsynonymous  mutation  or  a  premature  stop  codon  are  also  found  in  the  COSMIC  cancer  gene  census.  All  three  were 
found  in  either  high-grade  foci  or  metastatic  disease,  but  not  in  low-grade  foci.  GERP,  genomic  evolutionary  rate  profiling  score;  LN,  lymph  node 
metastasis;  NS,  nonsynonymous;  STOP,  introduction  of  premature  stop  codon;  *,  stop  codon. 


For  the  cases  in  which  there  was  evidence  of  a  common 
progenitor  for  both  low-  and  high-grade  disease  (PrCa  25  and 
PrCa  6),  the  nature  of  the  common  progenitor  is  not  clear. 
Histological  low-grade  cancer,  HGPIN  or  even  histologically 
uninvolved  prostate  are  all  possible  progenitors.  If  the  progenitor 
is  not  histological  low-grade  cancer,  high-grade  and  low-grade 
disease  may  have  emerged  simultaneously.  Characterization  of 
the  common  progenitor  requires  further  study. 

The  present  study  demonstrates  there  is  late  divergence  of 
high-grade  cancer  and  synchronous  metastatic  disease,  confirm¬ 
ing  high-grade  disease  as  the  precursor  to  metastasis  and  sug¬ 
gesting  few  if  any  additional  mutations  are  required  for 
malignant  cells  to  gain  the  ability  to  metastasize.  Given  the  dif¬ 
ficulty  in  obtaining  tissue  from  typical  bone  metastases  in  pros¬ 
tate  cancer,  more  easily  accessible  prostate  lesions  might  offer 
a  substitute  for  identifying  genetic  derangements  leading  to  pri¬ 
mary  resistance  of  metastatic  disease.  However,  mechanisms  of 
secondary  resistance  are  likely  induced  or  enriched  by  the 
selective  pressure  of  therapy  and  thus  are  likely  found  at  low 
levels,  if  at  all,  in  the  primary  lesion. 

The  most  common  gene  fusion  event  in  prostate  cancer  is  the 
TMPRSS2-ERG  fusion,  leading  to  high  levels  of  expression  of 
ERG,  which  is  otherwise  expressed  at  very  low  levels  in  prostate 
cancer.  We  evaluated  for  the  presence  of  this  fusion  through 
IHC  for  ERG  and  those  findings  supported  the  evidence  from 
the  exome  sequencing  data.  A  limitation  of  exome  sequencing  is 
the  inability  to  identify  additional  chromosomal  rearrangements. 
Although  this  is  a  limitation  it  was  not  required  for  the  goal  of 
the  present  study,  which  was  to  identify  tumor  lineages. 

One  goal  of  sequencing  prostate  tumors  is  to  identify  bio¬ 
markers  that  will  supplement  the  prognosis  given  by  clinical 
features  alone.  To  that  end  we  sought  to  identify  mutated 
genes  or  pathways  that  were  enriched  in  high  grade  and  meta¬ 
static  disease.  TP53  is  the  most  frequently  mutated  gene  in 
advanced  prostate  cancer  and  we  found  the  p53  signaling 
pathway  to  harbor  mutations  only  in  high  grade  and  metastatic 
disease.  Thus,  members  of  the  p53  signaling  pathway  are 


intriguing  candidates  for  potential  biomarkers  of  aggressive 
disease.  Recently,  expression  of  a  panel  of  cell  cycle  progres¬ 
sion-associated  genes  has  been  correlated  with  outcomes  of 
those  with  low-  or  intermediate-grade  prostate  cancer.(32)  We 
also  found  enrichment  for  cell  cycle-related  genes  among  those 
mutated  in  high-grade  disease,  although  this  was  of  borderline 
statistical  significance  in  this  small  cohort. 

Given  the  good  outcomes  of  low-grade  prostate  cancer,  we 
hypothesized  that  mutations  not  found  in  low-grade  disease  are 
more  likely  to  be  drivers  of  progression  and  metastasis  and 
therefore  more  important  potential  therapeutic  targets.  Indeed, 
the  three  genes  identified  here  that  are  also  found  in  the  COS¬ 
MIC  cancer  gene  census  were  identified  exclusively  in  high 
grade  and  metastatic  foci.  Of  note,  all  three  mutations  are  in 
highly  conserved  regions  as  evidenced  by  the  high  genomic 
evolutionary  rate  profiling  score  and  thus  are  likely  to  be 
biologically  significant.  Identification  of  driver  genes  through 
this  strategy  has  the  potential  to  guide  the  development  of 
future  therapies. 

The  present  study  contributes  to  the  growing  body  of  work 
demonstrating  that  next-generation  sequencing  from  microdis- 
sected  FFPE  prostatectomy  specimens  is  feasible.  Based  on 
these  techniques,  we  found  that  coincident  low-  and  high-grade 
prostate  cancer  can  emerge  through  either  independent 
progression  or  early  divergence  from  a  common  progenitor. 


Fig.  5.  Model  of  prostate  cancer  evolution.  Low-grade  and  high- 
grade  cancer  foci  progress  largely  in  parallel,  diverging  early  from  a 
common  progenitor.  Metastatic  disease  exhibits  late  divergence  from 
high-grade  disease.  It  cannot  be  ruled  out  that  in  some  cases  low-  and 
high-grade  diseases  arise  through  independent  origins. 
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High-grade  disease  is  characterized  by  mutations  in  known 
cancer-associated  genes  and  the  p53  signaling  pathway.  Meta¬ 
static  disease  is  closely  related  to  high-grade  primary  foci. 
These  data  can  be  leveraged  to  help  identify  potential  biomar¬ 
kers  and  drivers  of  progression  of  clinically  significant  disease. 
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Next-gen  tissue:  preservation  of  molecular  and 
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Abstract:  Background:  Personalization  of  cancer  therapy  requires  molecular  evaluation  of  tumor  tissue.  Traditional 
tissue  preservation  involves  formalin  fixation,  which  degrades  the  quality  of  nucleic  acids.  Strategies  to  bank  frozen 
prostate  tissue  can  interfere  with  diagnostic  studies.  PAXgene  is  an  alternative  fixative  that  preserves  protein  and 
nucleic  acid  quality.  Methods:  Portions  of  prostates  obtained  from  autopsy  specimens  were  fixed  in  either  10% 
buffered  formalin  or  PAXgene,  and  processed  and  embedded  in  paraffin.  Additional  sections  were  immediately 
embedded  in  OCT  and  frozen.  DNA  and  RNA  were  extracted  from  the  formalin-fixed,  PAXgene-fixed,  or  frozen  tissue. 
Quantitative  PCR  was  used  to  compare  the  quality  of  DNA  and  RNA  obtained  from  all  three  tissue  types.  In  addition, 
5  pm  sections  were  cut  from  specimens  devoid  of  cancer  and  from  prostate  cancer  specimens  obtained  at  prosta¬ 
tectomy  and  fixed  in  PAXgene.  They  were  either  stained  with  hematoxylin  and  eosin  or  interrogated  with  antibodies 
for  p63,  PSA  and  p504.  Results:  Comparable  tissue  morphology  was  observed  in  both  the  formalin  and  PAXgene- 
fixed  specimens.  Similarly,  immunohistochemical  expression  of  the  P63,  PSA  and  P504  proteins  was  comparable 
between  formalin  and  PAXgene  fixation  techniques.  DNA  from  the  PAXgene-fixed  tissue  was  of  similar  quality  to  that 
from  frozen  tissue.  RNA  was  also  amplified  with  up  to  8-fold  greater  efficiency  in  the  PAXgene  fixed  tissue  compared 
to  the  formalin-fixed  tissue.  Conclusions:  Prostate  specimens  fixed  with  PAXgene  have  preserved  histologic  morphol¬ 
ogy,  stain  appropriately,  and  have  preserved  quality  of  nucleic  acids.  PAXgene  fixation  facilitates  the  use  of  prosta¬ 
tectomy  tissue  for  molecular  biology  techniques  such  as  next-generation  sequencing. 

Keywords:  Prostate  cancer,  PAXgene,  formalin,  fixation,  next-generation  sequencing,  immunohistochemistry 


Introduction 

Prostate  cancer  (PCa)  is  the  most  common 
malignancy  in  many  Western  countries,  affect¬ 
ing  over  230,000  men  each  year  in  the  U.S. 
alone,  and  it  is  the  second  leading  cause  of 
cancer  death  among  American  men  [1]. 
Localized  PCa  is  a  clinically  heterogeneous  dis¬ 
ease  with  significant  variation  in  patient  out¬ 
come.  Determining  the  most  beneficial  man¬ 
agement  strategy  requires  risk  stratification 
based  on  clinical  factors  and  diagnostic/prog¬ 
nostic  biomarkers  [2].  Though  advanced  PCa  is 
also  clinically  heterogeneous,  most  men  are 
treated  with  the  same  therapies.  Biomarker- 
driven  therapies  are  just  now  beginning  to  enter 
into  clinical  trials,  and  none  have  entered  into 
routine  clinical  practice. 

It  is  essential  to  identify  better  predictors  that 
could  guide  therapeutic  actions  and  help  reach 


the  target  of  personalized  medicine  in  PCa.  In 
this  challenging  context,  next-generation  sequ¬ 
encing  (NGS)  technologies  constitute  powerful 
tools  to  identify  potential  prognostic  and  predic¬ 
tive  biomarkers  based  on  each  patient’s  muta¬ 
tional  landscape  and  expression  signature.  The 
recent  rise  of  sequencing  of  PCa  genomes  has 
significantly  increased  our  understanding  of  the 
molecular  basis  of  the  disease  [3-6].  Efforts  are 
being  made  towards  translating  these  new  find¬ 
ings  into  clinical  care  [7,  8]. 

The  development  of  biomarkers  through  NGS  or 
other  methods  requires  access  to  high  quality 
specimens.  In  the  case  of  PCa,  this  resource  is 
particularly  difficult  to  acquire.  PCa  is  rarely  vis¬ 
ible  on  gross  examination  due  to  its  small  size 
or  the  intermingling  of  benign  and  malignant 
glands,  making  sampling  difficult  and  often 
poorly  representative  of  the  tumor.  Moreover, 
standard  prostatectomy  specimen  processing 
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Figure  1.  Autopsy  prostate  processing.  Mirrored  sec¬ 
tions  from  four  prostate  specimens  from  autopsies 
were  fixed  either  in  PAXgene  or  formalin.  Two  speci¬ 
mens  were  fixed  in  PAXgene  for  24  hours,  and  two 
specimens  were  fixed  in  PAXgene  for  48  hours.  For 
each  prostate  a  frozen  sample  was  taken  before  fixa¬ 
tion. 


requires  that  the  gland  should  be  fully  fixed  in 
formaldehyde  solution  before  sectioning  to 
guarantee  the  accurate  assessment  of  surgical 
margins  [9].  Others  have  proposed  methods  to 
freeze  a  larger  fraction  of  the  specimen,  but 
these  methods  compromise  the  “capsule”  [10]. 
Formalin  fixation  causes  the  formation  of  chem¬ 
ical  cross  links  of  nucleic  acids  and  proteins, 
thus  resulting  in  degraded,  fragmented  DNA 
and  RNA  [11].  The  current  study  sought  to  eval¬ 
uate  an  alternative,  non-crosslinking  fixation 
reagent  for  prostate  specimens.  We  demon¬ 
strate  that  fixation  of  prostates  with  PAXgene 
results  in  preservation  of  high  quality  DNA  and 
RNA  with  preserved  tissue  architecture  and 
staining  characteristics. 

Materials  and  methods 

Tissue  preparation 

Four  random  prostates  devoid  of  PCa  were 
obtained  at  autopsy  and  prepared  by  the 
Human  Tissue  Resource  Center  (HTRC)  at  the 
University  of  Chicago.  Before  fixation,  a  punch 
biopsy  was  immediately  embedded  in  OCT  com¬ 
pound  and  frozen.  The  remainder  of  each  speci¬ 
men  was  bisected.  Half  was  fixed  in  10%  buff¬ 
ered  formalin  for  24  hours  using  standard 
fixation  techniques  and  embedded  in  paraffin. 
The  other  half  was  placed  in  PAXgene  Tissue 
FIX,  150  ml_,  or  adjusted  to  ~4  ml_  per  gram  of 
tissue,  for  either  24  or  48  hours  (Figure  1). 
Specimens  were  then  incubated  in  ethanol- 
diluted  PAXgene  Tissue  STABILIZER,  150  mL, 
or  adjusted  to  ~4  mL  per  gram  of  tissue,  for  at 
least  three  hours.  PAXgene-fixed  specimens 


were  embedded  in  paraffin  using  standard  eth¬ 
anol-based  techniques,  avoiding  all  exposure  to 
formalin. 

Similarly,  as  part  of  an  Internal  Review  Board- 
approved  protocol,  specimens  with  PCa  from 
four  patients  undergoing  radical  prostatectomy 
were  fixed  in  PAXgene  Tissue  FIX  for  48  hours 
and  in  PAXgene  Tissue  STABILIZER  for  at  least 
three  hours.  They  were  then  embedded  in  par¬ 
affin  using  standard  ethanol-based  tech¬ 
niques. 

Hematoxylin  and  eosin  (H&E)  staining  and  im- 
munohistochemistry  (IHC) 

Five  pm  sections  from  PAXgene-fixed  paraffin 
embedded  (PAXPE)  and  formalin-fixed  paraffin 
embedded  (FFPE)  sample  blocks  were  stained 
with  hematoxylin  and  eosin  according  to  stan¬ 
dard  procedures.  Two  pathologists  with  exper¬ 
tise  in  genitourinary  malignancies  (T.A.,  G.P.) 
performed  independent  morphology  assess¬ 
ment  including  overall  morphology  and  nuclear, 
cytoplasmic  and  cell  membrane  details  of  the 
prostatic  gland  and  stroma,  while  blinded  to  the 
fixation  method. 

For  IHC,  tissue  sections  were  treated  and 
stained  using  the  automated  VentanaBench- 
Mark  XT  System  (PSA)  or  the  Leica  Bond  III 
System  (P63  and  P504).  For  PSA  detection, 
anti-PSA  antibody  (DAKO,  cat#A0562,  Rabbit 
IgG,  polyclonal,  1:1500)  was  applied  on  tissue 
sections  and  the  antigen-antibody  binding  was 
detected  with  Ultraview  detection  kit  from 
Ventana.  For  P63  detection,  tissue  sections 
were  treated  with  Bond  Epitope  Retrieval 
Solution  2  (AR9640,  Leica,  Biosystems,  EDTA, 
pH  9)  for  20  minutes.  Anti-P63  antibody 
(Biocaremedical,  cat#CM163,  clone:  4A4, 
mouse  lgG2a,  1:100)  was  applied  on  tissue 
sections  for  25  minutes.  For  P504  detection, 
tissue  sections  were  treated  with  Bond  Epitope 
Retrieval  Solution  1  (AR9961,  Leica,  Biosy¬ 
stems,  citrate  buffer,  pH  6)  for  30minutes.  Anti- 
P504S  antibody  (Biocaremedical,  cat#ACA- 
200BK,  Rabbit  IgG,  polyclonal,  1:50)  was 
applied  on  tissue  sections  for  25  minutes.  For 
P63  and  P504  the  antigen-antibody  binding 
was  detected  using  Leica  Bond  Refine  polymer 
detection  system  (Leica  Biosystems,  DS9800). 
Tissue  sections  were  briefly  immersed  in  hema¬ 
toxylin  for  counterstaining.  Images  were 


1228 


Am  J  Transl  Res  2015;7(7):1227-1235 


Preservation  of  prostate  tissue  in  fidelity 


Figure  2.  Preservation  of  morphology  and  antigenicity  depending  on  fixation  method.  H&E  staining  (left),  and  P63 
(center)  and  PSA  (right)  immunohistochemistry  of  prostate  tissues  from  autopsy  specimens  fixed  in  PAXgene  (top) 
for  48  hours  or  in  formalin  (bottom).  Magnification  X10. 


Figure  3.  Prostate  cancer  staining.  H&E  staining  (A)  and  P63  (B),  P504S  (C)  and  PSA  (D)  immunohistochemical 
staining  of  prostatectomy  tissues  fixed  in  PAXgene  for  48  hours  (top)  or  formalin  (bottom).  Magnification  X10. 


obtained  on  a  DM6000B  microscope  (Leica)  at 
lOx  magnification. 

Nucleic  acids  extraction 

Genomic  DNA  and  total  RNA  from  PAXgene- 
fixed  samples  were  isolated  using  PAXgene 
Tissue  DNA  kit  (PreAnalytiX)  and  PAXgene 
Tissue  miRNA  kit  (PreAnalytiX),  respectively. 
Genomic  DNA  and  total  RNA  from  formalin-fixed 
samples  were  isolated  using  QIAamp  DNA  FFPE 
Tissue  kit  (Qiagen)  and  miRNeasy  FFPE  kit 


(Qiagen),  respectively.  Genomic  DNA  and  total 
RNA  from  frozen  OCT  embedded  samples  were 
isolated  using  DNeasy  Blood  and  Tissue  kit 
(Qiagen)  and  miRNeasy  Mini  kit  (Qiagen), 
respectively.  All  extractions  were  performed 
using  the  manufacturer’s  recommended  proto¬ 
cols.  Extracted  DNA  concentration  was  mea¬ 
sured  by  Qbit®  fluorometric  quantification  (Life 
Technologies).  DNA  fragmentation  and  RNA 
integrity  number  (RIN)  were  evaluated  by 
Screentape®  electrophoresis  system  (Agilent 
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turer’s  instructions.  Reactions 
were  performed  on  a  ViiA  7 
Real-Time  PCR  System  (Life 
Technologies).  Thermal  cycling 
conditions  were  as  followed:  an 
activation  step  of  15  min  at 
95°C  and  40  cycles  of  15  s  at 
94°C,  30  s  at  62°C  and  30  s  at 
70°C.  PCR  reactions  were  run 
in  triplicate.  The  acceptance  cri¬ 
teria  for  single  reactions  were 
that  the  cycle  threshold  (CT) 
must  be  below  40  and  the  stan¬ 
dard  deviation  in  CT  between 
triplicate  reactions  must  be 
below  0.15.  Melting  curves  had 
to  be  free  of  extraneous  peaks. 
The  significance  of  the  differ¬ 
ence  between  PAXPE  and  FFPE 
samples  was  tested  using  a 
paired  Student’s  T-test. 


Figure  4.  DNA  fragment  length  depending  on  fixation  method.  Representa¬ 
tive  results  of  DNA  fragments  length  between  prostate  tissue  samples  fixed 
in  PAXgene  for  48  hours  or  in  formalin  evaluated  by  Screentape®  electro¬ 
phoresis  system.  The  electrophoregram  for  the  FFPE  sample  shows  smaller 
more  degraded  DNA  as  compared  to  the  PAXPE  sample. 


Results 

Preservation  of  tissue  morphol¬ 
ogy  and  immunohistochemistry 
characteristics 


Technologies).  After  extraction,  DNA  was  stored 
at  4°C  and  RNA  at-80°C. 

Reverse  transcription  assay 

For  each  sample,  1  pg  of  RNA  was  transcribed 
into  cDNA  using  Superscript  VILO  MasterMix 
(Life  Technologies)  according  to  the  manufac¬ 
turer’s  instructions.  After  reverse  transcription, 
cDNA  was  stored  at  -20°C. 

Real  time  PCR  assay 

For  testing  DNA  integrity,  one  forward  primer 
and  six  reverse  primers  (Supplementary  Table 
1)  were  designed  to  amplify  six  different  tran¬ 
script  sequences  with  lengths  of  88  to  817 
base-pairs  (bp)  of  the  GAPDH  gene.  Fortesting 
cDNA  integrity,  one  forward  primer  and  five 
reverse  primers  (Supplementary  Table  1)  were 
designed  to  amplify  five  different  transcript 
sequences  with  lengths  of  95  to  983  bp  of  the 
(3-Actin  (ACTB)  gene.  Real  time  quantitative 
PCR  assays  were  performed  with  10  ng  of 
genomic  DNA  and  one  microliter  of  1:4  cDNA 
dilution,  using  iTaq  Universal  SYBR  Green 
Supermix  (Bio-rad)  according  to  the  manufac- 


To  evaluate  the  preservation  of  tissue  morphol¬ 
ogy  under  PAXgene  fixation,  sections  of  forma¬ 
lin-fixed,  paraffin-embedded  (FFPE)  and  PAX- 
gene-fixed,  paraffin-embedded  (PAXPE)  speci¬ 
mens  were  stained  with  hematoxylin  and  eosin. 
On  pathology  assessment,  the  preservation  of 
tissue  morphology  in  specimens  was  similar 
regardless  of  fixation  agent  used  (Figure  2A).  It 
was  noted  that  using  equivalent  staining  proto¬ 
cols,  PAXgene-fixed  tissue  exhibits  stronger 
eosinophilic  stain,  which  does  not  affect  the 
nuclear  detail.  Some  compromise  in  morpholo¬ 
gy  was  noted  in  the  autopsy  sections  fixed 
using  either  technique,  which  was  likely  due  to 
extended  post-mortem  time  prior  to  fixation. 

Analysis  of  prostate  tissue  relies  on  immuno- 
histochemical  analysis  using  antibodies  to 
Prostate  Specific  Antigen  (PSA),  the  basal  cell 
marker  P63,  and  P504,  a  marker  supporting  a 
diagnosis  of  malignancy.  Using  standard  IHC 
workflows,  sections  of  PAXPE  and  FFPE  pros¬ 
tates  devoid  of  cancer  were  evaluated  for  PSA 
and  P63.  Comparable  staining  characteristics 
and  staining  intensity  were  observed  between 
PAXPE  and  FFPE  specimen  for  PSA  and  P63, 
(Figure  2B,  2C).  PSA  staining  was  detected  as  a 
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Figure  5.  Preservation  of  DNA  and  RNA  depending  on  fixation  method.  A. 
Preservation  of  DNA.  Line  plot  showing  the  difference  in  cycle  threshold  (Ct) 
in  real-time  PCR  targeting  GAPDH  gene  DNA,  between  PAXgene-fixed  and 
formalin-fixed  prostate  tissue  samples,  and  the  corresponding  Frozen  tissue. 
Each  line  represents  the  average  values  from  2  prostate  specimens  from 
autopsies  fixed  for  48  hours.  *:  statistical  significance,  P-value  <  0.05.  B. 
Preservation  of  RNA.  Line  plot  showing  threshold  cycle  (Ct)  of  real-time  RT- 
PCR  targeting  (3-Actin  gene  cDNAfrom  PAXgene-fixed  or  formalin-fixed  pros¬ 
tate  tissue  samples.  Each  line  represents  the  average  values  from  2  prostate 
specimens  from  autopsies  fixed  for  48  hours.  *:  statistical  significance,  P- 
value  <  0.05. 


very  intense  stain  in  luminal  cells  of  prostate 
glands  whereas  P63  showed  as  strong  nuclear 
staining  in  the  basal  cells. 

To  confirm  that  there  was  minimal  effect  of 
PAXgene  fixation  on  immunohistochemical 
analysis  of  tissue  containing  PCa,  four  prosta¬ 
tectomy  specimens  were  fixed  in  PAXgene 
reagent  for  48  hours.  As  seen  with  non-PCa  tis¬ 
sue,  H&E  staining  demonstrated  that  tissue 
morphology  and  architecture  was  well  pre¬ 
served  (Figure  3A).  PCa  specimens  underwent 
immunohistochemical  evaluation  for  P504,  in 
addition  to  PSA  and  P63.  All  three  stained 
strongly,  demonstrating  high  levels  of  immune 


reactivity  on  PAXPE  sections 

(Figure  3B-D). 

Preservation  of  nucleic  acid 
quality 

To  investigate  the  quality  of 
nucleic  acids  from  the  differ¬ 
ent  fixation  methods,  total 
DNA  and  RNA  were  extracted 
from  each  autopsy  prostate 
specimen.  Sample  yield  var¬ 
ied  highly  between  specimens 
but  did  not  show  significant 
differences  between  the  pres¬ 
ervation  methods  (data  not 
shown).  For  DNA  from  both 
types  of  tissue  the  average 
fragment  length  was  5Kb  or 
more.  PAXgene-fixed  samples 
were  -25%  longer  than  forma¬ 
lin-fixed  ones  (Figure  4).  For 
RNA  samples,  RIN  scores 
were  below  2,  and  there  was 
no  significant  difference  bet¬ 
ween  PAXgene-  and  Formalin- 
fixed  tissues.  Quality  of  the 
RNA  extracted  from  frozen  tis¬ 
sue  sample  from  each  autop¬ 
sy  specimen  ranged  from  RIN 
1.5  to  3. 

The  relative  quality  of  DNA 
and  RNA  was  evaluated  by 
real-time  PCR  and  reverse 
transcription  real-time  PCR, 
respectively.  To  measure  the 
DNA  template  quality,  as  a 
function  of  the  fixation  meth¬ 
od,  six  primer  pairs  were  designed  to  amplify 
fragments  of  increasing  length  of  the  GAPDH 
gene,  from  88  bp  to  817  bp.  Consistently,  DNA 
from  PAXPE  samples  demonstrated  4-8  fold 
increased  amplification  efficiency  (correspond¬ 
ing  to  2-3  lower  cycle  threshold  value)  com¬ 
pared  to  FFPE  samples  (Figure  5A). 

Five  primer  pairs  were  used  to  evaluate  cDNA 
quality,  producing  amplicons  of  increasing 
length  of  the  (3-Actin  house-keeping  gene,  from 
95  bp  to  983  bp.  While  RIN  scores  were  low  in 
both  formalin  and  PAXgene  tissues,  there  was 
a  significant  difference  in  cDNA  template  quali¬ 
ty,  with  up  to  8-fold  increased  amplification  effi¬ 
ciency  in  PAXPE  samples  (Figure  5B). 
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Figure  6.  Preservation  of  DNA  and  RNA  depending  on  PAXgene  fixation  time. 
A.  Preservation  of  DNA.  Line  plot  showing  the  difference  in  cycle  threshold 
(Ct)  in  real-time  PCR  targeting  GAPDH  gene  DNA,  between  prostate  tissue 
samples  fixed  in  PAXgene  for  24  or  48  hours,  and  the  corresponding  Frozen 
tissue.  Each  line  represents  the  average  values  from  2  prostate  specimens 
from  autopsies.  *:  statistical  significance,  P-value  <  0.05.  B.  Preservation 
of  RNA.  Line  plot  showing  threshold  cycle  (Ct)  of  real-time  RT-PCR  targeting 
p-Actin  gene  cDNA  from  prostate  tissue  samples  fixed  in  PAXgene  for  24  or 
48  hours.  Each  line  represents  the  average  values  from  2  prostate  speci¬ 
mens  from  autopsies.  *:  statistical  significance,  P-value  <  0.05. 


Fixation  with  PAXgene  for  48  hours  requires 
alterations  to  the  typical  tissue  processing 
workflow  and  a  delay  in  availability  of  tissue  for 
diagnostic  studies.  Therefore  we  also  evaluat¬ 
ed  the  quality  of  DNA  and  RNA  from  PAXPE 
prostate  specimens  fixed  for  24  hours.  As  seen 
in  Figure  6,  quality  of  both  DNA  and  RNA  from 
these  specimens  was  slightly  inferior  to  that 
fixed  for  48  hours. 

Discussion 

We  report  here  that  PAXgene  fixation  and  stabi¬ 
lization  reagents  preserve  morphology  and 
immunoreactivity  equivalent  to  formalin  fixa¬ 
tion,  but  with  superior  preservation  of  nucleic 
acid  integrity.  This  technique  can  facilitate  col¬ 
lection  of  tissue  with  high  quality  DNA  and  RNA, 
without  compromising  tissue  architecture  by 


disrupting  the  capsule  prior  to 
fixation,  and  without  limiting  tis¬ 
sue  collection  to  small  samples 
obtained  without  prior  full  path¬ 
ological  review.  This  represents 
a  particular  advantage  for  stud¬ 
ies  evaluating  the  entire  gland 
or  correlating  morphologic  or 
staining  features  with  molecu¬ 
lar  characteristics. 

Hematoxylin  and  eosin  stained 
PAXPE  prostate  tissues  provid¬ 
ed  good  representation  of  his¬ 
tologic  features,  similar  to  that 
of  formalin.  As  previously  des¬ 
cribed  [12,  13]  an  EC  FP7  proj¬ 
ect  aimed  to  improve  pre  ana¬ 
lytic  procedures,  the  PAXgene 
Tissue  System  (PAXgene)  both 
hematoxylin  and  eosin  were 
modestly  more  intense  in 
PAXPE  samples  than  in  the  cor¬ 
responding  FFPE  samples. 
Immunohistochemical  analysis 
revealed  that  fixation  with 
PAXgene  did  not  alter  immuno¬ 
reactivity  of  antigens  relevant 
for  PCa  as  compared  to  forma¬ 
lin  fixation.  Staining  and  coun- 
terstaining  were  intense  for  all 
markers  tested.  Importantly, 
immunohistochemical  analysis 
for  this  study  did  not  require 
specific  optimization  of  routine 
staining  procedures  for  PAXPE 
samples,  in  contrast  to  what  has  been  previ¬ 
ously  reported  [12]. 

Consistent  with  previous  studies,  we  found  that 
PAXgene  reagents  preserved  both  DNA  and 
RNA  quality  significantly  better  than  formalin 
[14-16].  Using  real-time  PCR  on  DNA  and  RNA 
obtained  from  prostate  autopsy  specimens,  we 
showed  that  PAXPE  samples  have  lower  CT  val¬ 
ues  compared  to  FFPE  samples.  This  difference 
increased  with  amplicon  length,  indicating  that 
PAXgene  allows  efficient  amplification  even  of 
long  fragments.  Chromosomal  rearrangements 
and  fusion  genes  appear  to  play  an  important 
role  in  prostate  cancer  development  and  pro¬ 
gression  [4,  17-19].  Identification  of  these 
structural  variants  using  next-generation 
sequencing  is  improved  by  sequencing  across 
longer  DNA  inserts,  thus  increasing  the  physi- 
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cal  coverage  of  the  genome  [20].  Sequencing 
libraries  made  from  FFPE  require  short  inserts 
to  allow  for  amplification  across  cross-linked 
DNA.  DNA  from  PAXPE  tissue  facilitates  use  of 
sequencing  libraries  more  amenable  to  identi¬ 
fying  the  rearrangements  characteristic  of 
prostate  cancer. 

We  demonstrate  here  significant  differences  in 
downstream  assays  using  RNA  extracted  from 
PAXPE  versus  FFPE  prostatectomy  samples, 
though  RNA  from  both  tissue  types  were  of  low 
quality.  It  should  be  noted  that  RNA  is  likely  to 
have  substantially  degraded  during  post-mor¬ 
tem  time  before  fixation,  which  is  supported  by 
the  low  quality  of  RNA  from  matched  frozen 
tissue. 

A  limitation  of  this  study  is  that  a  portion  of  the 
study  used  autopsy  specimens,  which  have 
already  undergone  some  degree  of  post-mor¬ 
tem  degradation  prior  to  PAXgene  and  formalin 
fixation.  This  strategy  was  undertaken  to  avoid 
compromising  the  capsule  of  clinical  samples 
prior  to  fixation.  Prostatectomy  specimens  fixed 
wholly  in  PAXgene  also  demonstrated  pre¬ 
served  morphology  and  staining  characteris¬ 
tics.  Although  we  demonstrate  superior  preser¬ 
vation  of  nucleic  acids  with  PAXgene,  the 
question  of  how  long  those  features  remain 
stable  is  still  to  be  answered. 

The  prolonged  duration  of  fixation  (48  hours) 
may  be  problematic  for  the  workflow  of  some 
centers.  Though  24  hour  fixation  is  inferior  to 
48  hour  fixation  in  preserving  DNA  quality,  it  is 
superior  to  formalin  and  preserves  staining 
characteristics  and  therefore  is  a  suitable  alter¬ 
native.  In  addition,  the  routine  use  of  PAXgene 
fixation  increases  financial  costs.  If  prohibitive¬ 
ly  expensive  for  routine  use,  PAX  gene  fixation 
is  an  excellent  method  for  preserving  speci¬ 
mens  of  high  research  interest. 

In  this  study  we  thus  demonstrate  the  novel 
use  of  PAXgene  for  fixation  of  prostatectomy 
specimens.  Its  use  improves  the  quality  of 
molecular  analysis  without  compromising  histo- 
pathological  analysis.  PAXgene  fixation  may 
prove  especially  useful  for  assays  for  which  lon¬ 
ger  amplicons  are  required,  such  as  for  chro¬ 
mosomal  rearrangements,  or  for  studies  mak¬ 
ing  use  of  more  tissue  than  can  be  obtained 
with  a  simple  punch. 
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Conclusions 

Given  the  difficulty  of  identifying  PCa  on  gross 
examination  and  the  potential  for  compromise 
of  the  prostate  “capsule"  if  tissue  is  divided 
prior  to  fixation,  banking  significant  amounts  of 
PCa  tissue  samples  is  challenging.  The  alterna¬ 
tive,  non-crosslinking  fixation  agent  PAXgene 
preserves  the  quality  of  DNA  and  RNA  isolated 
from  fixed  tissue  without  compromising  the 
quality  of  tissue  histology  or  staining  character¬ 
istics.  Incorporation  of  PAXgene  fixation  into 
banking  protocols  should  facilitate  the  wider 
availability  of  high  quality  PCa  tissues  for  next- 
generation  sequencing  and  other  studies. 
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Supplementary  Table  1.  Sequences  of  primers  used  for  real  time  PCR  assay 


Name 

Sequence  (5’-3’) 

GAPDH-F 

AGTCCCCAGAAACAGGAGGT 

GAPDH-R-88  bp 

GAPDH-R-270  bp 

GAPDH-R-420  bp 

GAPDH-R-486  bp 

GAPDH-R-693  bp 

GAPDH-R-817  bp 

B-ACTIN-F 

AGAGCGCGAAAGGAAAGAA 

GTTAGTCACCG  G  CAG  G  CTTT 

C  AC  AAG  AG  G  ACCTCCATAAACC 

ACCCATGACTCAGCTTCTCC 

GTTT  CCG  G  AAG  ACG  G  AAT  G 

G  G  G  AG  CACAG  GTAAGT  G  CAT 

AAAT  CTG  G  CACCACACCTT  C 

B-ACTIN-R-95  bp 

B-ACTIN-R-626  bp 

B-ACTIN-R-703  bp 

B-ACTIN-R-828  bp 

B-ACTIN-R-983  bp 

GTTG  G  CCTT  G  G  G  GTT  CAG 

AGAGGCGTACAGGGATAGCA 

AGGGCATACCCCTCGTAGAT 

TA  AT  GT  C  ACG  CACG  ATTT  CC 

AG  G  CTG  G  AAG  AGTG  CCT  CAG 
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Abstract 

Background:  Genetic  analysis  of  advanced  cancer  is  limited  by  availability  of  representative  tissue.  Biopsies  of 
prostate  cancer  metastasized  to  bone  are  invasive  with  low  quantity  of  tumor  tissue.  The  prostate  cancer  genome 
is  dynamic,  however,  with  temporal  heterogeneity  requiring  repeated  evaluation  as  the  disease  evolves.  Circulating 
tumor  cells  (CTCs)  offer  an  alternative/liquid  biopsy",  though  single  CTC  sequencing  efforts  are  laborious  with  high 
failure  rates. 

Methods:  We  performed  exome  sequencing  of  matched  treatment-naive  tumor  tissue,  castrate  resistant  tumor  tis¬ 
sue,  and  pooled  CTC  samples,  and  compared  mutations  identified  in  each. 

Results:  Thirty-seven  percent  of  CTC  mutations  were  private  to  CTCs,  one  mutation  was  shared  with  treatment-naive 
disease  alone,  and  62%  of  mutations  were  shared  with  castrate-resistant  disease,  either  alone  or  with  treatment-naive 
disease.  An  acquired  nonsense  mutation  in  the  Retinoblastoma  gene,  which  is  associated  with  progression  to  small 
cell  cancer,  was  identified  in  castrate  resistant  and  CTC  samples,  but  not  treatment-naive  disease. This  timecourse  cor¬ 
related  with  the  tumor  acquiring  neuroendocrine  features  and  a  change  to  neuroendocrine-specific  therapy. 

Conclusions:  These  data  support  the  use  of  pooled  CTCs  to  facilitate  the  genetic  analysis  of  late  stage  prostate 
cancer. 

Keywords:  Circulating  tumor  cells,  Castrate  resistant  prostate  cancer, Tumor  evolution,  Neuroendocrine  prostate 
cancer 


Background 

Prostate  cancer  biology  is  initially  dominated  by  activity 
of  the  androgen  receptor  (AR),  and  androgen  depriva¬ 
tion  therapy  (ADT)  is  the  backbone  of  therapy  for  those 
with  metastatic  prostate  cancer.  Castrate  resistance  is 
nearly  universal,  however,  requiring  additional  manage¬ 
ment  decisions.  The  heterogeneity  of  primary  prostate 
cancer  is  well  established  [1-3],  and  polyclonal  seeding 
of  metastatic  sites,  and  metastasis-to-metastasis  spread, 
appear  to  be  possible,  if  not  common  [4].  Moreover,  the 
prostate  cancer  genome  is  dynamic,  with  different  clones 
dominating  over  time  in  response  to  different  lines  of 
therapy  [5,  6].  As  the  number  of  approved  agents  that 
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prolong  survival  increases,  the  dynamic  clonal  nature 
of  this  disease  is  likely  to  become  a  greater  issue.  For 
example,  there  appears  to  be  a  new  entity,  intermediate 
between  conventional  adenocarcinoma  and  classic  neu¬ 
roendocrine  disease,  that  arises  with  resistance  to  newer 
AR-targeted  agents  [7]. 

Circulating  tumor  cells  (CTCs)  and  cell  free  DNA 
(cfDNA)  are  emerging  as  non-invasive  alternative  means 
to  interrogate  the  genetics  of  late  stage  disease  rather 
than  invasive  and/or  risky  biopsy  of  metastatic  disease. 
They  are  complimentary  methods,  each  with  advantages 
and  disadvantages.  The  dilute  tumor  fraction  of  cfDNA 
in  most  patients  requires  very  high  sequencing  cover¬ 
age  to  detect  mutations,  typically  limiting  analysis  to  a 
small  number  of  targeted  regions  [5].  On  the  other  hand, 
amplifying  DNA  from  single  circulating  tumor  cells  is 
inconsistent  and  low  yield.  Successful  efforts  from  indi¬ 
vidual  cases  have  sequenced  four  of  99  collected  CTCs, 
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or  pooled  the  results  from  19  separately  sequenced 
CTCs  [8,  9].  We  present  here  whole  exome  sequencing  of 
pooled  CTCs,  with  matched  treatment- naive  and  castrate 
resistant  tissue,  demonstrating  identification  in  pooled 
CTCs  of  clinically  relevant  mutations  acquired  late  in  the 
course  of  disease. 

Methods 

Patient  recruitment  and  patient  samples 

The  patient  provided  written  informed  consent  and  was 
enrolled  in  the  “Prostate  Cancer  Sample  Collection”  pro¬ 
tocol  of  the  University  of  Chicago  Medical  Center,  which 
was  approved  by  the  University  of  Chicago  ethics  com¬ 
mittee  (approval  reference  number  13-1295).  Biopsy 
specimens  were  processed  per  clinical  protocols  includ¬ 
ing  hematoxylin  and  eosin  stain  and  immunohistochemi- 
cal  analysis  of  PSA  expression.  Circulating  tumor  cells 
were  collected  as  previously  described  [10].  In  short, 
mononuclear  cells  were  enriched  from  15  ml  peripheral 
blood  and  stained  with  an  Alexa-488  conjugated  EpCAM 
antibody  (Biolegend,  1:100)  and  a  QDot800  conjugated 
CD45  antibody  (Invitrogen,  1:100).  CTCs  were  isolated 
by  FACS-sorting  EpCAM+/CD45—  cells  on  a  MoFlo 
XDP  flow-sorting  machine.  One  thousand  WBCs  were 
isolated  by  FACS-sorting  EpCAM— /CD45+  cells  as  rep¬ 
resentative  of  germline  DNA. 

Whole  genome  amplification  and  exome  sequencing 

Isolated  EpCAM+/CD45—  CTCs  were  divided  into  three 
equal  samples  of  500  cells  each  and  independently  sub¬ 
jected  to  whole  genome  amplification  (WGA)  through 
multidisplacement  amplification  using  REPEI-G  (Qia- 
gen),  following  the  manufacturer  s  instructions,  as  were 
EpCAM— /CD45+  WBCs  representative  of  germline. 
The  quality  of  amplification  was  evaluated  through  PCR 
amplification  of  eight  targets  [11]  to  evaluate  uniform 
amplification  across  multiple  chromosomes,  and  length 
of  amplification  product.  Primers  used  are  listed  in  Addi¬ 
tional  file  1.  Exome  sequencing  libraries  were  prepared 
using  NEB  Next  Ultra  (New  England  Biolabs)  kit  follow¬ 
ing  the  manufacturers  instructions.  Sequencing  was  per¬ 
formed  on  a  HiSeq  2000  (Illumina). 

Sequencing  analysis 

Illumina  sequencing  data  was  mapped  to  the  GRCh37 
version  of  the  human  reference  genome  using  BWA- 
MEM  [12]  and  further  processed  following  the  GATK 
Best  Practices  [13].  Somatic  variants  were  called  for  each 
tumor  or  CTC  sample  in  a  paired  manner  using  MuTect 
[14]  in  high  confidence  mode.  Potential  germline  vari¬ 
ants  were  removed  by  excluding  positions  with  germline 
coverage  <20  x,  variant  reads  in  the  germline,  or  variants 
present  in  dbSNP,  and  remaining  variants  were  enriched 


for  highest  confidence  by  including  only  variants  in  the 
exome  and  with  5  or  more  supporting  reads  comprising 
10%  or  more  of  all  reads.  Somatic  mutations  passing  all 
filters  (listed  in  Additional  file  2)  in  a  given  sample  were 
then  examined  for  the  presence/absence  of  supporting 
reads  in  other  samples  using  samtools  mpileup  [15,  16]. 
Mutations  were  visually  verified  by  examining  supporting 
reads  using  Alview  [17],  and  all  somatic  mutations  were 
annotated  using  AVIA  [18]. 

Results 

A  70  year  old  patient  was  diagnosed  with  metastatic  pros¬ 
tate  cancer,  confirmed  with  biopsy  of  a  bone  metastasis 
(Fig.  la).  He  underwent  androgen  deprivation  therapy 
but  eventually  developed  castrate  resistant  prostate  can¬ 
cer  (CRPC).  Additional  conventional  therapies  yielded 
short-lived  responses  and  subsequent  progression  of  dis¬ 
ease  (Fig.  lb).  A  liver  biopsy  performed  22  months  after 
the  initial  diagnosis  confirmed  metastatic  prostate  cancer 
(Fig.  lc),  but  with  neuroendocrine  and  small  cell  features, 
which  was  not  appreciated  on  the  initial  biopsy.  Given 
the  change  in  dominant  histology  he  started  treatment 
with  carboplatin  and  etoposide,  but  after  two  cycles  he 
entered  hospice  care. 

To  evaluate  the  extent  to  which  the  current  disease 
biology  can  be  discovered  through  CTCs,  CTCs  were  col¬ 
lected  from  15  ml  of  peripheral  blood  within  a  month  of 
the  liver  biopsy.  CTCs  were  enriched  by  FACS-sorting 
EpCAM+/CD45—  cells  as  previously  described  [10]  and 
divided  into  three  equal  fractions.  Multidisplacement 
amplification  was  used  to  amplify  the  entire  genome 
of  each  fraction  independently.  To  evaluate  for  loss  of 
coverage  due  to  lack  of  amplification  and  for  length  of 
amplified  fragments,  the  amplified  CTC  genomes  were 
evaluated  with  endpoint  PCR  for  eight  targets  of  vary¬ 
ing  lengths  across  six  different  chromosomes  (Fig.  Id). 
More  targets  were  able  to  be  amplified  from  one  CTC 
pool  (CTC1)  than  the  other  two  pools  (CTC2,  CTC3).  All 
three  samples  were  used  for  exome  sequencing. 

Whole  exome  sequencing  was  performed  on  a  ger¬ 
mline  sample  and  five  tumor  samples:  initial  diagnos¬ 
tic,  treatment-naive  biopsy  of  an  ischial  metastasis;  liver 
biopsy  of  CRPC;  and  three  samples  of  pooled  CTCs 
divided  from  a  single  CTC  collection.  The  fractions  of 
PCR  duplicates  were  all  0.25  or  less  (Fig.  2a).  Adequate 
coverage  was  obtained  for  the  germline  sample,  the  tis¬ 
sue  samples,  and  CTC1,  samples  that  were  predicted  to 
perform  well  by  initial  evaluation,  with  coverage  of  over 
98%  of  the  genome,  median  depths  of  44-106  x,  and 
10  x  coverage  of  90%  or  more  of  the  exome  (Fig.  2b).  The 
lower  quality  CTC  samples  covered  just  37  and  40%  of 
the  exome  at  10  x.  Due  to  their  poor  quality  they  were 
left  out  of  downstream  analysis.  Mutation  patterns  were 
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Fig.  1  Clinical  timecourse  of  disease,  a  Hematoxylin-Eosin  stain  (fop)  and  immunohistochemical  analysis  of  PSA  protein  ( bottom )  of  treatment-naive 
bone  biopsy  diagnosing  metastatic  prostate  cancer,  b  Clinical  timecourse  indicating  PSA  levels,  timing  of  biopsies,  timing  of  the  single  CTC  collec¬ 
tion  ~3  weeks  after  CRPC  biopsy,  which  was  divided  into  three  pools,  and  timing  of  therapies.  CAB  combined  androgen  blockade,  abi  abiraterone, 
doc  docetaxe I,  RT  palliative  radiation  therapy,  CE  carboplatin  and  etoposide.  c  CT  image  of  liver  metastasis  that  underwent  biopsy  {white  arrowhead). 
d  Fraction  of  PCR  targets  that  were  amplified  from  each  whole  genome  amplification  product 


consistent  in  the  treatment-naive,  CRPC,  and  CTC  sam¬ 
ple,  including  nucleotide  substitutions  (Fig.  2c)  and  type 
of  mutation  (Fig.  2d). 

The  allele  frequency  was  comparable  for  mutations  in 
the  tissue  specimens,  but  higher  in  the  CTCs  (Fig.  3a). 
Moreover,  while  trunk  mutations  (those  found  in  both 
treatment-naive  and  CRPC  tissue  samples)  were  near 
50%  allele  frequency,  suggesting  close  to  100%  pure  CTCs 
all  with  one  allele  mutated,  branch  (shared  with  one  tis¬ 
sue  sample)  and  leaf  (found  only  in  CTCs)  mutations 
were  at  lower  allele  frequency,  indicating  genetic  hetero¬ 
geneity  among  CTCs  (Fig.  3b). 

There  was  considerable  genetic  heterogeneity  among 
tumor  tissue  samples,  with  just  8%  of  mutations  being 
trunk  mutations,  and  25%  being  found  in  two  samples. 
Thirty-six  percent  of  mutations  from  CRPC  were  not 
initially  identified  in  the  treatment-naive  sample  but  had 
at  least  two  reads  supporting  their  presence,  indicating 


enrichment  in  the  CRPC  tissue  of  subclonal  populations 
from  the  treatment-naive  tissue.  The  reverse— mutations 
identified  in  the  treatment-naive  sample  with  low  fre¬ 
quency  supporting  reads  in  CRPC — was  not  identified 
(Fig.  3c,  d).  The  pooled  CTCs  identified  71%  of  muta¬ 
tions  shared  by  treatment-naive  and  CRPC  tissue  sam¬ 
ples,  suggesting  most  mutations  at  high  allele  frequency 
(early  mutations)  were  represented  in  the  CTCs  (Fig.  3e). 
Examining  the  two  low  quality  CTC  pools,  each  pool 
identified  22%  of  trunk  mutations,  consistent  with  their 
lOx  coverage  being  37-40%  of  the  genome.  All  the  trunk 
mutations  identified  in  the  low  quality  pools  were  also 
identified  in  the  high  quality  pool,  suggesting  consistency 
in  the  mutations  able  to  be  identified  by  CTCs. 

The  CTCs  had  the  highest  fraction  of  mutations  shared 
with  at  least  one  other  sample,  at  64%  (Fig.  3f).  Fifty  per¬ 
cent  were  trunk  mutations,  shared  with  both  tissue  sam¬ 
ples,  and  an  additional  14%  were  shared  with  the  CRPC 
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sample  alone.  A  single  mutation  was  shared  between 
CTCs  and  treatment-naive  tissue  and  not  found  in  CRPC 
tissue.  Though  the  other  two  CTC  samples  had  poor 
exome  coverage  and  were  not  of  sufficient  quality  to 
identify  mutations  a  priori,  33%  of  mutations  identified  in 
the  high  quality  CTC  sample  were  also  identified  in  one 
of  the  low  quality  CTC  samples. 

One  of  the  mutations  shared  between  the  CRPC  tis¬ 
sue  sample  and  CTCs  was  a  premature  stop  codon  in  the 
Retinoblastoma  (RBI)  gene  (Fig.  4a),  which  is  associated 
with  isolated  bilateral  retinoblastoma  and  meningioma 
[19].  The  appearance  of  this  RBI  mutation  coincided  with 
loss  of  adenocarcinoma  features,  including  PSA  expres¬ 
sion  (compare  Fig.  4b,  top  and  middle  panels,  to  Fig.  la), 
and  gain  of  neuroendocrine  features  in  the  tumor,  which 
was  confirmed  by  synaptophysin  expression  (Fig.  4b, 
bottom  panel).  This  phenotypic  change  toward  therapy- 
emergent  neuroendocrine  prostate  cancer  prompted  a 
change  in  management  strategy  away  from  standard  of 
care  for  CRPC. 

Discussion 

The  prostate  cancer  genome  is  heterogeneous,  both 
between  and  within  the  multiple  foci  characteristic  of 
primary  disease.  The  clonal  architecture  of  advanced 


disease  is  dynamic,  with  new  clones  gaining  dominance 
in  response  to  new  therapies.  Combined,  these  neces¬ 
sitate  repeated  genomic  or  molecular  assessments  over 
the  course  of  disease  in  order  to  have  a  complete,  current 
understanding  of  a  patient  s  personalized  disease. 

Circulating  tumor  cells  offer  a  non-invasive  mecha¬ 
nism  to  repeatedly  evaluate  the  shifting  dynamics  of 
disease.  This  has  demonstrated  clinically  meaningful 
evaluations  of  specific  alterations  [20,  21].  However, 
there  have  been  few  direct  correlations  between  treat¬ 
ment-naive  tissue,  CTCs,  and  contemporaneous  CRPC 
tissue,  and  wider  scale  genomic  evaluations  are  still  in 
early  stages.  Lohr  et  al.  used  whole  genome  sequencing 
to  evaluate  quality  of  amplified  DNA  followed  by  exome 
sequencing  of  19  individual  CTCs  to  demonstrate  late 
divergence  of  CTCs  and  a  previously  resected  lymph 
node  [8].  Jiang  et  al.  used  laser  capture  microdissection 
to  capture  and  evaluate  99  individual  CTCs  collected 
over  five  blood  collections,  identifying  four  individual 
CTCs  with  high  quality  DNA  and  eight  with  moder¬ 
ate  quality  [9].  The  CTCs  in  that  study  identified  15%  of 
trunk  SNVs,  with  supporting  reads  for  an  additional  14% 
of  reads. 

Using  our  pooled  CTC  strategy,  we  generated  suc¬ 
cessful  sequencing  libraries  from  33%  of  samples.  The 
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Fig.  3  Circulating  tumor  cells  pool  captures  majority  of  trunk  mutations  and  additional  mutations  from  metastatic  disease,  a  Variant  allele  frequen¬ 
cies  for  treatment-nai've,  CRPC,  and  CTC1  mutations,  b  Average  variant  allele  frequencies  of  trunk,  branch,  and  leaf  mutations  in  two  tissue  samples 
and  one  CTC  sample,  c  Allele  fraction  in  CRPC  tissue  of  mutations  identified  in  treatment-nai've  tissue,  and  d  vice  versa.  On  the  left  are  mutations 
identified  in  both  samples.  On  the  right  are  mutations  identified  in  only  sample.  In  d,  many  of  the  mutations  not  identified  independently  in  treat¬ 
ment-naive  tissue  had  evidence  they  were  present  at  low  allele  frequency,  e  Number  of  trunk  mutations  identified  in  three  pooled  CTC  sequencing 
libraries,  f  Fraction  of  CTC  mutations  shared  with  other  tissue  samples 


sequencing  demonstrated  high  correlation  with  tis¬ 
sue  samples,  confirming  the  biological  relevance  of 
the  CTC  exomes,  identifying  71%  of  trunk  mutations, 
along  with  additional  mutations  acquired  later  in  dis¬ 
ease.  This  includes  a  clinically  meaningful  mutation 
in  RBI  which  likely  contributed  to  a  change  in  pheno¬ 
type  to  neuroendocrine  features,  prompting  a  change 
in  management  strategies.  The  RBI  gene  is  altered  in 
nearly  9%  of  advanced  prostate  cancer  cases,  through 
deletion,  frameshift  mutations,  and  introductions  of 
premature  stop  codons  [22].  Beltran  et  al.  compared 
advanced  prostate  neuroendocrine  and  adenocarcinoma, 


demonstrating  that  RBI  alterations  are  significantly 
enriched  in  advanced  prostate  cancer  with  neuroendo¬ 
crine  features  (70%  altered)  compared  to  that  with  pure 
adenocarcinoma  features  (32%  altered)  [23].  Loss  of  RBI 
function  is  common  in  primary  small  cell  cancer  of  the 
prostate  or  lung,  and  in  animal  models  it  promotes  devel¬ 
opment  of  small  cell  carcinoma  [24,  25]. 

The  clonal  relationship  among  all  three  specimens  sug¬ 
gest  that  neuroendocrine  disease  arose  from  adenocarci¬ 
noma,  rather  than  being  a  coincident,  independent  clone. 
In  addition,  the  high  frequency  of  mutations  in  CRPC  tis¬ 
sue  that  were  present  at  low  frequency  in  treatment-nai’ve 
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Fig.  4  Identification  of  RBI  nonsense  mutation  in  CRPC  tissue  and  CTC  samples,  a  Integrative  Genomics  Viewer  [28]  images  of  the  S816X  mutation 
identified  in  CRPC  and  CTC1,  but  not  treatment-nai've  sample.  The  brown  bars  in  the  middle  panels  represent  the  mutated  nucleotide  and  its  posi¬ 
tion  within  the  sequencing  reads.  The  histograms  in  the  upper  panels  summarize  the  fraction  of  reference  {gray  or  blue )  or  mutated  {brown)  reads,  b 
Hematoxylin-Eosin  stain  {top)  and  immunohistochemical  analysis  of  PSA  {middle)  or  synaptophysin  protein  {bottom)  of  the  CRPC  liver  biopsy 


tissue  supports  the  idea  that  advanced  disease,  includ¬ 
ing  neuroendocrine  disease,  arises  from  subclonal 
population(s)  in  the  initial  specimen.  Of  note,  though 
CTCs  from  patients  with  neuroendocrine  prostate  cancer 
are  more  frequently  nonclassical  than  those  with  patients 
with  adenocarcinoma  (17),  the  RBI  mutation  was  identi¬ 
fied  in  classical  EpCAM-f-  CTCs. 

Fewer  mutations  were  identified  in  CTCs  than  in  treat- 
ment-nai’ve  or  CRPC  tissue  samples.  The  significance  of 
this  is  unclear.  It  may  be  that  the  limited  number  of  CTCs 
was  unable  to  capture  the  extensive  diversity  of  clones 
comprising  disease  in  the  tissue.  Alternatively,  it  may  be 
that  CTCs  represent  a  limited  number  of  aggressive  and 
clinically  relevant  clones.  The  CTCs  were  not  clonal,  as 
evidenced  by  the  presence  of  branch  and  leaf  mutations. 
This  genetic  heterogeneity  among  CTCs  is  supported  by 


Massard  et  al.  [26]  based  on  a  single  genomic  alteration, 
the  ERG  alteration  pattern.  The  extent  to  which  CTCs 
represent  all  the  relevant  subclones  needs  to  be  explored 
further. 

There  are  several  advantages  of  a  pooled  CTC  strat¬ 
egy  over  single  CTC  sequencing,  including  availability 
of  resources.  Our  strategy  relied  on  FACS -sorting  and 
whole  genome  amplification  using  a  commercially  avail¬ 
able  kit,  which  are  readily  available  to  most  research¬ 
ers.  We  did  not  require  laser  capture  microdissection  or 
robotic  micromanipulation.  Disadvantages  include  appli¬ 
cability  limited  to  patients  with  a  higher  burden  of  CTCs 
and  inability  to  fully  characterize  heterogeneity  at  the  cel¬ 
lular  level. 

While  we  had  a  much  higher  success  rate  sequencing 
pooled  CTCs  than  has  been  reported  with  single  CTC 
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sequencing,  only  one  of  three  pools  provided  high  qual¬ 
ity  data.  This  may  have  been  due  in  part  to  our  multiple 
displacement  amplification  strategy  [27].  This  was  chosen 
for  the  low  error  rate  of  its  polymerase,  but  it  may  be  less 
reliable  in  amplifying  the  majority  of  the  genome  com¬ 
pared  to  PCR-based  methods. 

As  the  number  of  effective  therapies  used  for  treat¬ 
ment  of  advanced  prostate  cancer  increases,  there  is 
an  increasing  appreciation  of  the  dynamic  nature  of 
the  genomics  of  advanced  disease  in  response  to  thera¬ 
peutic  pressure.  We  demonstrate  here  that  sequencing 
pooled  CTCs  is  a  feasible,  noninvasive,  and  informa¬ 
tive  way  to  evaluate  the  current  molecular  features  of 
advanced  disease. 

Conclusions 

The  histology,  behavior,  and  genomics  of  advanced  pros¬ 
tate  cancer  evolve  in  response  to  therapeutic  pressure. 
Pooled  CTCs  are  a  feasible,  non-invasive  way  to  interro¬ 
gate  the  molecular  characteristics  of  advanced  prostate 
cancer. 
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Dear  Dr.  Ganz, 

We  would  like  to  submit  the  manuscript  "Genomic  heterogeneity  within  individual 
prostate  cancer  foci  impacts  predictive  biomarkers  of  targeted  therapy"  for 
consideration  for  publication  as  an  Article  in  JNCI. 

Most  lethal  prostate  cancer  progresses  from  relapse  of  primary  disease.  Given  known 
heterogeneity  between  multiple  foci,  prognosis  is  generally  attributed  to  a  dominant 
index  focus.  The  heterogeneity  within  the  index  focus  of  aggressive,  non-metastatic 
disease,  however,  is  not  known. 

We  use  multiregion  genomic  analysis  to  examine  the  heterogeneity  of  potentially  lethal, 
treatment-naive  prostate  cancer.  Our  data  demonstrate  several  significant  and 
clinically  important  findings.  First,  there  is  marked  genomic  heterogeneity  in  high  risk, 
treatment-naive  disease,  indicating  the  heterogeneity  found  in  metastatic  castrate 
resistant  prostate  cancer  (mCRPC)  can  also  be  found  at  much  earlier  stages.  Second, 
this  heterogeneity  extends  to  potential  genetic  biomarkers,  such  as  defects  in 
homologous  recombination  repair  or  alterations  in  the  PI3K  pathway.  Biomarkers 
present  in  one  region  are  often  discordant  with  other  regions,  and  biomarkers  identified 
in  the  prostate  tumor  are  discordant  with  synchronous  metastases.  Third,  there  are 
numerous  co-occurring  alterations  in  pathways  typically  associated  with  mCRPC, 
though  found  in  subclonal  populations  in  treatment-naive  disease.  Fourth,  there  is 
widespread  dispersal  of  subclonal  populations,  so  even  widespread  alterations  are 
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often  subclonal. 

Currently  there  is  uncertainty  about  how  to  manage  high  risk,  non-metastatic  disease, 
in  part  due  to  uncertainty  about  its  biology.  The  LATITUDE  and  STAMPEDE  trials  are 
practice  changing  to  include  abiraterone  with  prednisone  for  castrate  sensitive  disease 
with  extensive  metastases.  Physicians  are  less  comfortable  extending  this 
management  to  high-risk,  non-metastatic  disease,  though  these  patients  were  included 
in  the  STAMPEDE  trial.  Though  there  is  clearly  a  failure-free  survival  benefit  for  non¬ 
metastatic  disease,  the  overall  survival  data  is  not  mature.  Moreover,  there  is 
uncertainty  regarding  the  biology  of  nonmetastatic  disease  relative  to  metastatic 
disease. 

The  data  in  this  study  are  the  most  thorough  genomic  characterization  to  date  of 
potentially  lethal,  nonmetastatic,  treatment-naive  disease.  While  the  findings 
presented  here  cause  concern  about  the  potential  for  resistance  to  systemic  therapy 
that  is  associated  with  such  marked  genetic  heterogeneity,  they  also  demonstrate  that 
the  biology  of  aggressive,  non-metastatic  disease  is  similar  to  metastatic,  therapy- 
resistant  disease.  Given  that  biological  similarity,  these  data  give  the  rationale  for  the 
patient  population  of  the  STAMPEDE  trial,  pooling  high  risk  localized  disease  (as 
included  in  the  present  study)  with  widely  metastatic  disease.  And  they  lend  support  to 
embrace  treatment  of  non-metastatic  disease  in  a  similar  fashion  to  widely  metastatic 
disease. 

Though  PARP  inhibitors  and  Akt  inhibitors  are  not  yet  approved  for  prostate  cancer, 
they  are  the  targeted  therapies  on  the  horizon  for  mCRPC.  One  response  to  the 
CHAARTED,  STAMPEDE,  and  LATITUDE  trials  is  the  desire  to  move  additional 
systemic  therapies  to  an  earlier  disease  state.  The  data  presented  here  also  provide 
insight  into  using  biomarkers  from  intraprostatic  disease  to  guide  use  of  those 
therapies. 

Our  manuscript  has  significant  implications  for  the  tens  of  thousands  of  patients 
diagnosed  with  potentially  lethal,  aggressive  primary  prostate  cancer  every  year.  We 
feel  this  timely  manuscript  will  be  welcomed  by  the  prostate  cancer  community  as  it 
deliberates  about  how  best  to  manage  high  risk,  treatment-naive  disease.  We  hope 
you  agree. 

Thank  you  for  your  consideration  of  this  work  for  publication  in  JNCI.  Please  feel  free 
to  contact  me  directly  should  you  have  any  questions. 

Kind  regards, 

David  J.  VanderWeele  MD,  PhD 

Assistant  Clinical  Investigator 

Laboratory  for  Genitourinary  Cancer  Pathogenesis 

LGCP,  NCI,  NIH 

37  Convent  Drive,  Rm  1066A 

Bethesda,  MD  20892 

301.443.2622 

david.vanderweele@nih.gov 
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Dear  Dr.  Ganz, 

We  would  like  to  submit  the  manuscript  “Genomic  heterogeneity  within  individual  prostate  cancer  foci 
impacts  predictive  biomarkers  of  targeted  therapy”  for  consideration  for  publication  as  an  Article  in  JNCI. 

Most  lethal  prostate  cancer  progresses  from  relapse  of  primary  disease.  Given  known  heterogeneity 
between  multiple  foci,  prognosis  is  generally  attributed  to  a  dominant  index  focus.  The  heterogeneity 
within  the  index  focus  of  aggressive,  non-metastatic  disease,  however,  is  not  known. 

We  use  multiregion  genomic  analysis  to  examine  the  heterogeneity  of  potentially  lethal,  treatment- 
naive  prostate  cancer.  Our  data  demonstrate  several  significant  and  clinically  important  findings. 

First,  there  is  marked  genomic  heterogeneity  in  high  risk,  treatment-naive  disease,  indicating  the 
heterogeneity  found  in  metastatic  castrate  resistant  prostate  cancer  (mCRPC)  can  also  be  found  at 
much  earlier  stages.  Second,  this  heterogeneity  extends  to  potential  genetic  biomarkers,  such  as 
defects  in  homologous  recombination  repair  or  alterations  in  the  PI3K  pathway.  Biomarkers  present  in 
one  region  are  often  discordant  with  other  regions,  and  biomarkers  identified  in  the  prostate  tumor  are 
discordant  with  synchronous  metastases.  Third,  there  are  numerous  co-occurring  alterations  in 
pathways  typically  associated  with  mCRPC,  though  found  in  subclonal  populations  in  treatment-naive 
disease.  Fourth,  there  is  widespread  dispersal  of  subclonal  populations,  so  even  widespread  alterations 
are  often  subclonal. 

Currently  there  is  uncertainty  about  how  to  manage  high  risk,  non-metastatic  disease,  in  part  due  to 
uncertainty  about  its  biology.  The  LATITUDE  and  STAMPEDE  trials  are  practice  changing  to 
include  abiraterone  with  prednisone  for  castrate  sensitive  disease  with  extensive  metastases. 

Physicians  are  less  comfortable  extending  this  management  to  high-risk,  non-metastatic  disease, 
though  these  patients  were  included  in  the  STAMPEDE  trial.  Though  there  is  clearly  a  failure-free 
survival  benefit  for  non-metastatic  disease,  the  overall  survival  data  is  not  mature.  Moreover,  there  is 
uncertainty  regarding  the  biology  of  nonmetastatic  disease  relative  to  metastatic  disease. 

The  data  in  this  study  are  the  most  thorough  genomic  characterization  to  date  of  potentially  lethal, 
nonmetastatic,  treatment-naive  disease.  While  the  findings  presented  here  cause  concern  about  the 
potential  for  resistance  to  systemic  therapy  that  is  associated  with  such  marked  genetic  heterogeneity, 
they  also  demonstrate  that  the  biology  of  aggressive,  non-metastatic  disease  is  similar  to  metastatic, 
therapy-resistant  disease.  Given  that  biological  similarity,  these  data  give  the  rationale  for  the  patient 
population  of  the  STAMPEDE  trial,  pooling  high  risk  localized  disease  (as  included  in  the  present 
study)  with  widely  metastatic  disease.  And  they  lend  support  to  embrace  treatment  of  non-metastatic 
disease  in  a  similar  fashion  to  widely  metastatic  disease. 


Though  PARP  inhibitors  and  Akt  inhibitors  are  not  yet  approved  for  prostate  cancer,  they  are  the 
targeted  therapies  on  the  horizon  for  mCRPC.  One  response  to  the  CHAARTED,  STAMPEDE,  and 
LATITUDE  trials  is  the  desire  to  move  additional  systemic  therapies  to  an  earlier  disease  state.  The 
data  presented  here  also  provide  insight  into  using  biomarkers  from  intraprostatic  disease  to  guide  use 
of  those  therapies. 

Our  manuscript  has  significant  implications  for  the  tens  of  thousands  of  patients  diagnosed  with 
potentially  lethal,  aggressive  primary  prostate  cancer  every  year.  We  feel  this  timely  manuscript  will 
be  welcomed  by  the  prostate  cancer  community  as  it  deliberates  about  how  best  to  manage  high  risk, 
treatment-naive  disease.  We  hope  you  agree. 

Thank  you  for  your  consideration  of  this  work  for  publication  in  JNCI.  Please  feel  free  to  contact 
me  directly  should  you  have  any  questions. 


Kind  regards, 


David  J.  YanderWeele  MD,  PhD 

Assistant  Clinical  Investigator 

Laboratory  for  Genitourinary  Cancer  Pathogenesis 

LGCP,  NCI,  NIH 

37  Convent  Drive,  Rm  1066 A 

Bethesda,  MD  20892 

240.760.6829 

david.  vanderweele  @  nih.gov 
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ABSTRACT 


Background:  Most  lethal  prostate  cancer  progresses  from  relapse  of  aggressive  primary  disease. 
Given  known  heterogeneity  between  multiple  foci,  prognosis  is  generally  attributed  to  a 
dominant  index  focus.  Targeted  therapies  show  promise  for  aggressive  prostate  cancer,  and 
predicted  response  is  based  on  the  presence  of  genomic  biomarkers.  The  present  study 
characterizes  the  genomic  heterogeneity,  including  that  of  predictive  biomarkers,  within  the 
index  focus  of  aggressive  primary  prostate  cancer. 

Methods:  Genomic  characterization  was  performed  on  70  regions  from  the  index  foci  of  ten 
patients  with  treatment-naive,  aggressive  primary  prostate  cancer.  Five  to  nine  regions  of  each 
index  focus  were  evaluated  by  exome  and  low  coverage  whole  genome  sequencing.  Intrafocal 
genomic  heterogeneity  and  heterogeneity  of  alterations  that  predict  response  to  therapy  was 
determined. 

Results:  Exome  sequencing  and  copy  number  estimates  demonstrate  branched  evolution  with 
>75%  of  point  mutations  being  subclonal,  including  numerous  pathways  associated  with  castrate 
resistant  prostate  cancer.  There  is  wide  dispersal  of  comingled  subclones,  with  individual 
subclonal  populations  detectable  across  all  regions  sampled.  Intraindividual  heterogeneous 
alterations  in  5  genes  predict  response  to  targeted  therapies.  Within  biomarker-positive  cases, 
25%  of  intraprostatic  regions  are  biomarker-negative,  with  discordance  between  intraprostatic 
regions  and  lymph  node  metastases. 

Conclusions:  Treatment-naive,  nonmetastatic  prostate  cancer  has  marked  intrafocal 
heterogeneity.  Numerous  alterations  in  pathways  associated  with  castrate  resistant  prostate 
cancer  are  present  in  subclonal  populations,  including  biomarkers  predictive  of  response  to 


targeted  therapy. 
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INTRODUCTION 

Thirty  percent  of  patients  receiving  definitive  therapy  for  newly  diagnosed  prostate 
cancer  relapse,  with  higher  rates  of  relapse  for  those  with  high  risk  or  locally  advanced  disease  1_ 
3.  For  over  seven  decades,  systemic  therapy  has  focused  on  the  androgen  receptor  (AR),  a 
disease-specific  target.  One  goal  of  genomic  characterizations  of  both  untreated  and  advanced 
prostate  cancer  is  to  identify  subtypes  of  disease  that  predict  sensitivity  to  new  therapies  4~6. 
These  studies  have  suggested  that  in  advanced  disease  over  60%  of  patients  have  targetable 
altered  pathways  in  addition  to  the  AR  signaling  pathway  6.  The  PARP  inhibitor  olaparib  is  the 
targeted  therapy  furthest  in  clinical  development,  for  patients  with  tumors  harboring  mutations  in 
DNA  repair  pathways,  comprising  approximately  20%  of  advanced  disease  6_8.  Recently 
ipatasertib  has  also  demonstrated  efficacy  in  patients  with  loss  of  the  tumor  suppressor  gene 
PTEN  9.  It  is  hoped  that  further  efforts  will  define  additional  subsets  of  disease  with  specific 
susceptibilities,  effectively  dividing  the  diagnosis  into  several  different  molecularly-based 
diagnoses. 

The  efficacy  and  durability  of  response  to  targeted  therapies  depends  in  part  on  the 
prevalence  of  the  targeted  alteration  within  the  tumor  cell  population.  Multifocal  prostate  cancer 
is  known  to  have  early  divergence  between  foci,  but  there  is  often  a  dominant,  high  grade  index 
focus  considered  to  represent  clinically  significant  disease  10_14.  Supporting  this,  there  is 
relatively  little  intrapatient  heterogeneity  in  metastatic  disease  when  prostate  cancer  becomes 
lethal 15,16.  However,  the  extent  of  heterogeneity  within  the  index  focus  of  aggressive  prostate 
cancer  is  unknown. 
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We  undertook  multiregion  genomic  evaluation  of  an  individual  index  focus  from  ten 
patients  with  localized  (lymph  node  negative)  or  locally  advanced  (lymph  node  positive)  disease, 
representing  those  most  likely  to  relapse.  To  our  knowledge,  this  is  the  first  genomic  evaluation 
focused  on  potentially  lethal,  locally  advanced  prostate  cancer.  The  present  study  describes 
branched  genomic  evolution  characterized  by  marked  heterogeneity  in  both  localized  and  locally 
advanced  disease.  Compared  to  localized  disease,  locally  advanced  disease  is  characterized  by 
increased  copy  number  changes  and  numerous  subclonal  alterations  in  pathways  associated  with 
castrate  resistant  prostate  cancer  (CRPC).  Significantly,  alterations  predictive  of  response  to 
targeted  therapy  are  heterogeneous  within  the  prostate  and  between  intraprostatic  disease  and 
lymph  node  metastatic  disease. 


METHODS 
Subjects  and  consent 

Prostatectomy  specimens  were  collected  from  subjects  enrolled  in  “Prostate  Cancer 
Sample  Collection”  tissue  collection  protocol  of  the  University  of  Chicago  Medical  Center  (see 
Supplemental  Methods).  All  individuals  included  in  the  study  provided  written  informed 
consent.  All  cases  were  reviewed  by  a  genitourinary  pathologist  (GP),  and  3D  reconstructions  of 
the  index  focus  were  created  based  on  location  and  histologic  appearance  of  tumor  within  each 
block. 

Genomic  analysis 

Germline  DNA  was  extracted  from  blood  using  a  DNeasy  kit  (Qiagen),  except  for  case 
1027,  when  it  was  obtained  from  histologically  benign  prostate  tissue.  DNA  was  extracted  from 


5 


6 


tumor  regions  using  PAXgene  Tissue  DNA  kit  (Qiagen)  and  exome  and  low  coverage  whole 
genome  sequencing  performed  (see  Supplemental  Methods).  A  subset  of  variants  was  verified 
with  ultradeep  targeted  resequencing  on  a  MiSeq  (Supplemental  Table  1).  Subsequently 
seventeen  regions  were  excluded  from  further  analysis  due  to  low  tumor  cell  purity  (<30%) 
based  on  variant  allele  frequency  of  identified  mutations  leaving  70  regions  with  sufficient 
sequencing  quality  and  tumor  cell  purity.  Positions  harboring  mutations  that  passed  quality  filters 
were  then  evaluated  for  all  regions  in  the  corresponding  case  using  samtools  mpileup.  Mutations 
were  considered  present  if  the  variant  allele  frequency  (VAF)  >1%,  and  VAF  >  10%  in  at  least 
one  region.  Indels  were  confirmed  using  IGV. 

Copy  number  measures  were  generated  by  BICseq2  17  (see  Supplemental  Methods) 
setting  filters  to  p-value  <  0.01  and  log2  ratio  <  -0.2  or  >  0.2  17 .  Whole  genome  sequencing  data 
was  unavailable  for  case  1009,  and  copy  number  alterations  were  determined  using  CNVkit 18 
where  log2  ratio  was  <  -0.2  or  >  0.2  and  confidence  interval  excluded  one.  To  avoid  falsely 
elevating  the  fraction  of  subclonal  events  due  to  lower  tumor  purity  in  some  samples,  if  log2 
ratios  in  all  regions  from  the  case  were  in  the  same  direction  (>0  or  <0),  an  alteration  was 
considered  clonal.  If  not,  it  was  considered  subclonal. 

Immunohistochemistry  was  performed  for  ERG  expression.  Nuclear  ERG  expression 
intensity  and  staining  distribution  were  scored  on  a  scale  of  0-3  by  a  genitourinary  pathologist 
(GP). 

3D  representations  of  prostatectomy  specimens  and  network  trees 
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Representations  of  prostatectomy  specimens  were  created  using  Povray  based  on  manual 
mapping  of  the  tumor  focus  and  regions  sampled.  Network  analysis  was  constructed  using 
Network  19,20  from  features  comprising  all  somatic  mutations  identified. 

Mutations  per  cell 

The  average  number  of  mutation  per  cell  was  calculated  from  exome  sequencing  data 
based  on  Landau  et  al  21  (see  Supplemental  Methods).  Average  number  of  mutations  per  cell 
was  calculated  as  (Mutation  VAF)  *  (purity  *  ploidy  +  2  *  (1  -  purity))  /  purity. 

Predictive  genetic  biomarkers 

Potential  biomarkers  were  queried  using  The  Drug  Gene  Interaction  Database 
(dgidb.genome.wustl.edu)  22.  A  biomarker  was  considered  positive  if  it  contained  both  a 
nonsynonymous  mutation  and  deletion  (log2  <  -0.2)  or  a  deep  deletion  (log2  <  -0.4).  To  avoid 
falsely  elevating  heterogeneity  of  biomarkers  a  deletion  of  log2  ratio  <  -0.2  was  also  considered 
a  biomarker  without  mutation  if  at  least  one  region  in  the  case  had  deep  deletion  at  that  locus. 

Statistics 

A  two-sided  student’s  t-test  using  R  was  used  to  analyze  the  difference  between  the 
number  of  mutations  and  the  percent  genome  altered  in  localized  versus  locally  advanced  cases. 
A  log-rank  Mantel-Cox  test  was  performed  using  Prism  to  evaluate  the  difference  in  biochemical 
recurrence  between  cases  with  more  or  less  than  the  median  CRPC-associated  alterations. 


RESULTS 
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Intrafocal  heterogeneity  of  aggressive  primary  prostate  cancer 

To  examine  the  heterogeneity  of  treatment-naive,  aggressive  primary  prostate  cancer, 
multiregion  sequencing  was  performed  on  an  individual  index  focus  from  the  prostatectomy 
specimen  of  ten  patients  (Supplemental  Figure  1).  Clinical  characteristics  are  typical  of 
patients  with  aggressive  disease,  with  Gleason  scores  7-9  and  pathologic  tumor  stage  pT2-T3 
(Table  1).  Five  cases  are  localized  (lymph  node  negative),  five  are  locally  advanced  (lymph  node 
positive).  Three  patients  had  biochemical  recurrence  following  prostatectomy,  and  four 
underwent  additional  therapy  after  prostatectomy.  For  each  case,  exome  and  low  coverage 
whole  genome  sequencing  was  performed  on  five  to  nine  regions  of  a  single  contiguous  index 
focus.  Lymph  node  metastases  obtained  at  the  time  of  prostatectomy  and  seminal  vesicle  regions 
involved  through  direct  extension  were  included  for  two  cases  each  (1010  and  1022;  1015  and 
1022,  respectively). 

We  identified  887  unique  mutation  events  across  70  index  focus  regions  (88.7 
variants/case,  ~1.3  nonsilent  variants/Mb)  (Supplemental  Table  2),  with  a  higher  number  of 
mutations  in  locally  advanced  than  in  localized  cases  (108.6  vs  68.6,  respectively).  A  subset  of 
mutations  was  verified  with  ultra-high  coverage  (~40,000x)  targeted  resequencing.  Of  107 
mutations  assessed,  2  (2%)  were  not  confirmed,  and  none  was  found  in  germline  (Supplemental 
Table  1).  Allele  frequencies  measured  by  exome  sequencing  were  similar  to  those  measured  by 
targeted  resequencing  (R2=0.861,  Supplemental  Figure  2). 

Exome  sequencing  demonstrated  marked  heterogeneity  of  single  nucleotide  mutations 
and  indels  within  the  cancer  genome  of  the  index  focus  in  all  cases  examined  (Figure  1A).  In 
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each  focus  trunk  mutations  (mutations  identified  in  all  regions)  comprised  the  minority. 
Alterations  in  ERG,  SPOP,  FOXA1,  IDH1  are  among  the  drivers  proposed  to  define  prostate 
cancer  subtypes  in  75%  of  cases  5.  Five  of  10  cases  presented  here  could  be  assigned  one  of 
these  subtypes  (3  ERG  overexpression,  1  FOXA1  mutation,  1  SPOP  mutation)  (Supplemental 
Figure  3).  As  with  a  single  region  sequencing  approach,  subtypes  were  mutually  exclusive.  As 
expected,  these  alterations  occurred  early,  with  the  alteration  identified  in  every  region  sampled 
in  all  but  one  case  (ERG  in  case  1029). 

Stability  of  mutation  processes  over  time 

Several  studies  have  presented  lines  of  evidence  for  a  field  effect  in  prostate  cancer  13,23, 
suggesting  there  is  a  long  accumulation  of  mutations  over  the  life  of  the  patient.  Given  this,  we 
examined  the  types  of  alterations  contributing  to  early  (trunk),  late  (leaf),  and  middle  (branch) 
events.  In  contrast  to  the  60%  or  more  clonal  mutations  shown  for  10  other  solid  tumors  24 ,  fewer 
than  25%  of  variants  identified  were  trunk  variants.  Despite  the  numerically  increased  number 
of  mutations  in  locally  advanced  cancer,  there  was  no  increase  in  fraction  of  late  mutations  in 
these  cases  (Supplemental  Figure  4A).  Unlike  other  malignancies  24_26,  we  also  found  no 
evidence  for  changes  in  mutational  processes  over  time,  with  similar  patterns  of  transitions  vs 
transversions  and  patterns  of  non-synonymous,  synonymous,  premature  stop,  and  splice  site 
mutations  among  trunk,  branch,  and  leaf  mutations  (Supplemental  Figure  4B,C). 

Heterogeneity  among  structural  changes 

Chromosomal  events  such  as  TMPRSS2:ERG  fusion  appear  to  be  early  events  in  prostate 
cancer,  present  even  in  a  fraction  of  premalignant  lesions  27 ,28.  To  identify  the  regional  pattern  of 
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chromosomal  changes,  we  performed  copy  number  estimates  from  low-coverage  whole  genome 
sequencing  (Figure  IB,  Supplemental  Table  3).  There  were  many  canonical  copy  number 
changes  that  were  estimated,  including  high  frequency  of  loss  of  chromosomes  8p,  13,  and  16, 
focal  deletion  on  chromosome  10,  and  gain  of  chromosomes  8q  and  7.  Locally  advanced  cases 
had  higher  mean  percent  genome  altered  than  those  without  lymph  node  involvement,  which  is 
associated  with  a  worse  prognosis  29,30  (Figure  1C).  Interestingly,  two  of  the  four  cases 
harboring  regions  with  the  highest  percent  of  the  genome  altered  (>  0.2)  also  contained  regions 
with  low  copy  number  burden  (-0.05),  suggesting  there  is  focal  genome  instability. 

Correlation  between  histology  and  genetic  features 

We  constructed  network  trees  for  each  case  based  on  mutation  patterns.  These  were 
overlaid  on  3D  representations  of  the  prostatectomy  specimens  (Figure  2,  Supplemental  Figure 
5).  In  cases  with  histologic  heterogeneity,  it  appeared  the  relationship  between  regions  was  as 
likely  to  correlate  with  histologic  similarity  than  spatial  proximity.  In  case  1010  for  example, 
despite  their  proximity  to  regions  C36A  and  C46A,  respectively,  regions  C36B  and  C46B  share 
high  grade  features  and  late  divergence  of  their  cancer  genome,  supported  by  shared  mutations 
and  copy  number  changes.  Both  are  more  similar  histologically  and  diverge  later  from  C50A 
than  from  C36A  and  C46A. 

We  also  observed  early  divergence  of  spatially  proximal  regions  without  obvious 
correlations  with  histology  (early  divergence  of  region  B3A  in  case  1024,  Supplemental  Figure 
5).  In  most  cases,  however,  there  was  consistent  histologic  appearance  across  the  focus 
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(Supplemental  Figure  1),  and  spatially  proximal  regions  showed  late  divergence.  This  was 
most  clearly  seen  in  cases  1024  and  1034  (Supplemental  Figure  5). 

Subclonal  alterations  in  clinically  significant  pathways 

A  recent  report  of  the  genomic  landscape  of  castrate  resistant  prostate  cancer  (CRPC) 
demonstrated  alterations  in  statistically  or  clinically  significant  genes  from  several  pathways  that 
contribute  to  prostate  cancer  progression,  including  androgen  receptor,  cell  cycle,  PI3K, 
RAF/RAS,  WNT,  DNA  repair,  and  chromatin  modifier  pathways  6.  Many  of  these  genes  are 
found  to  be  altered  at  low  frequency  in  treatment-naive  disease 5.  We  hypothesized  these  genes 
might  be  altered  in  more  cases  of  treatment-naive  disease  than  previously  appreciated,  but  in 
subclonal  populations  below  the  level  of  detection  in  other  studies.  Indeed,  all  ten  cases  have 
alterations  in  at  least  two  of  41  CRPC -related  genes  (Figure  3A).  Moreover,  those  with  higher 
than  median  subclonal  alterations  in  CRPC  pathways  were  more  likely  to  have  recurrence  than 
those  below  the  median  (p-value  <  0.05,  log-rank  test)  (Figure  3B).  The  majority  of  the 
alterations  we  identified  were  chromosomal  events,  highlighting  the  frequency  and  importance  of 
chromosome  level  events  compared  to  single  nucleotide  mutations  31 . 

Given  that  many  of  the  CRPC  pathways  altered  are  potential  therapeutic  targets,  we 
evaluated  the  heterogeneity  in  genetic  biomarkers  that  predict  response  to  targeted  therapy. 

Using  concurrent  copy  loss  and  mutation  or  deep  deletion  we  identified  five  genes  that  predict 
susceptibility  to  two  therapies  with  demonstrated  efficacy  in  CRPC  8’9,  with  individual  cases 
harboring  biomarker  alterations  in  up  to  three  separate  genes  (Figure  3C).  One  case  (1009)  had 
germline  mutation  in  the  DNA-binding  domain  of  BRCA2  that  is  of  unknown  significance 
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(p.Ala2730Pro),  with  concurrent  somatic  copy  loss  of  BRCA2.  The  remaining  predictive 
biomarkers  were  heterogeneous  among  regions  sampled  such  that,  of  biomarker-positive  cases, 
25%  of  regions  were  biomarker-negative.  Moreover,  there  was  discordance  between 
intraprostatic  regions  and  lymph  node  metastases:  in  two  cases  with  lymph  node  metastases 
available,  only  one  of  three  biomarker  alterations  found  in  intraprostatic  regions  was  also  found 
in  the  lymph  node,  and  one  of  two  alterations  identified  in  a  lymph  node  was  found  in  the 
prostate  (Figure  3C). 

Spatial  dispersion  of  subclones 

For  a  targeted  therapy  to  be  effective,  the  targeted  alteration  ought  to  be  found  in  every 
cancer  cell.  We  assume  that  mutations  detected  in  every  region  of  a  focus  (trunk  mutations) 
occur  early.  It  is  not  clear,  however,  if  these  mutations  are  detected  in  every  region  because  they 
are  found  in  every  cell,  or  because  they  are  found  in  subclones  that  have  dispersed  widely.  We 
examined  the  variant  allele  frequencies  (VAFs)  of  mutations  in  coincident  regions  to  determine 
if  wide  spatial  dispersion  indicates  it  is  indeed  an  early  event.  As  expected,  trunk  mutations  are 
found  at  high  allele  frequency,  generally  correlating  with  tumor  cell  purity  (typically  -0.30, 
representing  -60%  tumor  cell  purity;  Figure  4A).  Leaf  mutations,  detected  in  a  single  region, 
have  low  allele  frequency  within  that  region.  Most  branch  mutations  also  have  low  allele 
frequency  within  a  region,  indicating  subclones  can  have  wide  (present  in  multiple  regions)  but 
sparse  (low  allele  frequency)  dispersion. 

To  evaluate  this  further,  we  generated  2D  density  plots  for  representative  spatially 
distinct  regions,  representing  VAF  normalized  for  tumor  cell  purity  and  ploidy  (akin  to  cancer 
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cell  fraction)  (Figure  4B,C).  This  demonstrates  that  branch  subclone  populations  have  highly 
variable  cancer  cell  fractions,  consistent  with  a  selection  mode  of  tumor  evolution  32.  Notably, 
even  mutations  identified  in  all  regions  (trunk  mutations)  can  have  low  allele  frequency  within  a 
region,  indicating  they  are  found  in  a  widely  dispersed  subclone  rather  than  being  truly 
ubiquitous. 


DISCUSSION 

We  present  here  the  largest  multiregion  genomic  evaluation  of  the  index  focus  of 
aggressive  primary  prostate  cancer.  This  study  demonstrates  the  heterogeneity  of  treatment- 
naive,  aggressive  disease,  with  fewer  than  25%  of  mutations  found  in  all  regions  of  the  focus. 
Coincident  foci  are  known  to  be  divergent,  leading  to  reliance  on  the  dominant  index  focus  for 
prognostication.  These  data  demonstrate  remarkable  heterogeneity  within  the  clinically  relevant 
index  focus. 

The  treatment  of  prostate  cancer  has  recently  entered  the  genomic  era,  with  PARP 
inhibitor  therapy  based  on  alterations  in  DNA  repair  pathways  6-8  demonstrating  response  in 
88%  of  biomarker  positive  patients,  and  an  Akt  inhibitor  has  also  demonstrated  benefit  in 
biomarker-positive  patients  9.  There  is  interest  in  moving  these  therapies  up  to  an  earlier  disease 
state,  with  trials  underway  testing  olaparib  in  the  neoadjuvant  or  biochemically  recurrent  space 
(NCT02324998,  NCT03047135),  where  biomarker  evaluation  is  presumably  based  on  the  index 
lesion  of  the  primary  tumor.  In  the  present  study  we  show  that  while  prostate  driver  alterations 
are  early  events,  predictive  biomarkers  such  as  ATM,  BRCA1,  and  BRCA2  can  be  altered  in  a 
small  fraction  of  cells  of  the  index  focus.  Moreover,  we  find  that  even  wide  spatial  dispersal  of 
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an  alteration  does  not  imply  it  is  ubiquitous.  In  these  cases,  one  would  expect  a  transient 
response  to  targeted  therapy.  This  is  in  line  with  studies  demonstrating  that  intratumoral 
heterogeneity  is  linked  with  poor  outcome  and  therapy  resistance  25’33-35. 

Currently  there  is  uncertainty  about  how  to  manage  high  risk,  non-metastatic  disease,  in 
part  due  to  uncertainty  about  its  biology.  The  LATITIUDE  and  CHAARTED  trials  36,37  indicate 
a  remarkable  survival  benefit  in  moving  life-prolonging  therapies  up  to  an  earlier  disease  state, 
and  the  STAMPEDE  trial  38,39  suggests  this  also  applies  to  non-metastatic  disease,  though  the 
survival  data  is  less  mature.  The  findings  presented  here  indicate  that  aggressive,  non¬ 
metastatic  disease,  especially  locally  advanced  disease,  is  remarkably  similar  to  metastatic, 
castrate-resistant  disease.  It  harbors  a  high  percent  of  the  genome  altered  by  copy  number 
alterations,  and  there  are  numerous  alterations  typically  associated  with  metastatic  CRPC. 

Our  data  expands  the  growing  body  of  work  demonstrating  heterogeneity  of  primary 
prostate  cancer.  Previous  work  demonstrated  interfocal 10,11,40  or  intra-  and  interfocal 
heterogeneity  12,13,41,  though  these  have  focused  on  lower  risk,  lower  stage  disease,  a  disease 
state  that  is  likely  to  be  cured  by  surgery  and  is  more  dissimilar  to  widely  metastatic  disease. 
Multiple  expression-based  assays  have  been  developed  to  refine  prognosis  of  localized  disease, 
and  Wei  et  al41  demonstrated  intra-  and  to  a  greater  extent  interfocal  heterogeneity  in  these 
assays.  Two  studies42,43  have  compared  primary  and  metastatic  disease  and  demonstrated  cases 
of  biomarker-negative  prostate  tumors  and  biomarker-positive  metastases  due  to  the 
accumulation  of  additional  alterations  in  metastatic  disease,  which  would  presumably  lead  to  a 
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false  negative  primary  tissue  biomarker.  The  present  study  now  demonstrates  a  positive  prostate 
tissue  biomarker  not  representative  of  all  tumor  cells,  predicted  to  lead  to  poor  response. 

These  data  demonstrate  there  can  be  marked  heterogeneity  among  regions  of  an 
individual  index  focus  in  aggressive  primary  prostate  cancer.  This  heterogeneity  also  applies  to 
genetic  biomarkers,  highlighting  the  risk  of  selecting  biomarker-driven  therapies  based  on 
primary  tissue.  The  concern  about  the  heterogeneity  of  CRPC  has  led  to  the  proposal  of  novel 
therapeutic  strategies  to  compensate  for  this  44 .  Our  data  suggest  concerns  about  tumor 
heterogeneity  also  apply  to  the  management  of  androgen  sensitive  disease. 
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Figure  legends 

Figure  1.  Intrafocal  heterogeneity  of  index  foci.  A.  Exome  sequencing  was  performed  on  five 
to  nine  regions  of  a  single  contiguous  index  focus  from  ten  cases  of  treatment-naive  prostate 
cancer.  Presence  (green)  or  absence  (white)  of  mutations  within  each  region  sampled  is  shown. 
Regions  are  represented  by  columns,  mutations  by  rows.  Five  locally  advanced,  node  positive 
(LN  pos)  cases  are  above,  five  localized  (LN  neg)  cases  are  below.  Previously  identified 
prostate  cancer  driver  genes  are  indicated  in  black.  Upper  panels,  the  level  of  the  radical 
prostatectomy  (RP)  specimen  from  which  the  region  was  sampled,  from  white  (apex)  to  black 
(base),  with  purple  (seminal  vesicles)  and  red  (lymph  node  metastasis)  indicating  extraprostatic 
regions.  Left  panels  specify  mutations  as  trunk  (found  in  all  regions),  branch  (at  least  two 
regions  and  not  all  regions),  or  leaf  (identified  in  one  region)  mutations.  B.  Log2  plots  of  copy 
number  estimates  based  on  low  coverage  whole  genome  sequencing,  except  case  1009  based  on 
exome  sequencing.  Left  panel,  the  level  of  the  prostatectomy  specimen  from  which  the  region 
was  sampled,  as  in  A.  Right  panel,  percent  genome  altered  by  copy  number  variation.  C.  Mean 
fraction  of  the  genome  altered  in  locally  advanced  (LN  pos)  and  localized  (LN  neg)  cases.  Error 
bars  represent  standard  deviation. 


Figure  2.  Correlation  between  genomic  and  histologic  features.  Network  analysis  of  case 
1010  based  on  shared  mutations  among  eight  regions  of  a  contiguous  index  focus  and  one  region 
from  a  lymph  node  metastasis  is  overlaid  on  a  cartoon  representation  of  the  prostatectomy 
specimen.  Corresponding  H&E  images  of  each  region  that  underwent  exome  sequencing  are 
depicted  on  the  left.  Met  =  lymph  node  metastasis. 
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Figure  3.  Heterogeneity  in  CRPC-related  alterations  and  predictive  biomarkers.  A. 

Alterations  across  all  cases  in  a  panel  of  clinically  relevant  genes  found  to  be  altered  in  advanced 
prostate  cancer  6.  B.  Kaplan-meier  analysis  of  biochemical  recurrence  between  cases  with 
greater  or  fewer  than  the  median  CRPC-associated  alterations  (6),  p-value  <0.05  (log-rank  test). 
C.  Fraction  of  intraprostatic  regions  within  a  case  that  harbor  deletion  and  mutation  or  deep 
deletion  in  five  predictive  biomarker  genes,  fraction  of  regions  represented  by  size  of  the  blue 
circle  (deletion)  or  “M”  (mutation).  Two  cases  (1010,  1022)  had  lymph  node  metastases 
evaluated,  represented  by  gray  shading  of  the  square. 


Figure  4.  Dispersion  of  subclones.  A.  Kernal  density  plot  of  variant  allele  frequency  (VAF)  of 
trunk  (brown),  branch  (orange),  and  leaf  (yellow)  mutations  across  all  index  focus  mutations. 

B,C.  Calculated  average  number  of  mutations  per  diploid  genome,  or  cancer  cell  fraction,  for 
each  mutation  identified  in  representative  pairs  of  coincident  regions.  Each  point  represents  one 
mutation  (B,  case  1009;  C,  case  1010).  Colors  are  as  in  A.  D.  YAF,  relative  to  the  mean  VAF  of 
trunk  mutations,  of  early  prostate  cancer  driver  gene  mutations  for  all  regions  of  the  case.  Gray, 
diploid.  Blue,  one  copy  deleted. 
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Table  1 

Clinical  information  for  cases  evaluated 


1009 

1010 

1015 

1022 

1029 

1003 

1020 

1024 

1027 

1034 

Stage  category 

Loc  adv 

Loc  adv 

Loc  adv 

Loc  adv 

Loc  adv 

Localized 

Localized 

Localized 

Localized 

Localized 

Age  at  surgery 

57 

66 

57 

61 

51 

56 

62 

60 

60 

64 

Ethnicity 

AA 

White 

AA 

White 

White 

White 

White 

White 

White 

AA 

PSA  (ng/mL) 

8.35 

12.76 

43.05 

NA 

7.96 

4.91 

4.30 

14.20 

9.20 

5.00 

Pathologic  stage 

T3bNl 

T3bNl 

T3bNl 

T3bNl 

T3aNl 

T2cN0 

T2c 

T3aN0 

T3bN0 

T2aN0 

Gleason  grade 

4+5 

4+5 

3+4  (5) 

4+5 

4+3 

4+5 

3+4 

3+4 

4+3  (5) 

4+3 

Tumor  volume  (%) 

40 

60 

65 

45 

25 

15 

25 

25 

25 

20 

Surgical  margin 

positive 

negative 

positive 

positive 

negative 

negative 

negative 

positive 

negative 

negative 

Nodes  examined 

22 

20 

18 

14 

15 

5 

4 

22 

22 

13 

Nodes  involved 

6 

1 

2 

4+5 

1 

0 

0 

0 

0 

0 

Germline  sample 

blood 

blood 

blood 

blood 

blood 

blood 

blood 

blood 

benign 

prostate 

blood 

Prostate  or  S V 
regions 

7 

8 

8 

6 

8 

5 

9 

8 

7 

6 

Lymph  node 
regions 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

Time  to  recurrence 
(mos) 

1 

2 

3 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

Additional  therapy 

Yes 

Yes 

Yes 

Yes 

No 

No 

No 

No 

No 

No 
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